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ABSTRACT 

We present new results from the SCUBA Local Universe Galaxy Survey (SLUGS), the first 
^ \ large systematic submillimetre survey of the local Universe. Since our initial survey of a 

sample of 104 //M5-selected galaxies we have now completed a survey of a sample of 81 
optically-selected galaxies, observed with the SCUBA camera on the James Clerk Maxwell 
Telescope. Since SCUBA is sensitive to the 90% of dust too cold to radiate significantly in 
the IRAS bands our new sample represents the first unbiased survey of dust in galaxies along 
iy-j ' the whole length of the Hubble sequence. 

, We find little change in the properties of dust in galaxies along the Hubble sequence, 

except a marginally significant trend for early-type galaxies to be less luminous submillimetre 
sources than late-types. We nevertheless detected 6 out of 11 elliptical galaxies, although 
some of the emission may possibly be synchrotron rather than dust emission. As in our earlier 
Q ■ work on IRAS galaxies we find that the IRAS and submillimetre fluxes are well-fitted by a 

\ two-component dust model with dust emissivity index /3=2. The major difference from our 

C/3 i earlier work is that we find the ratio of the mass of cold dust to the mass of warm dust is 

' much higher for our optically-selected galaxies and can reach values of ^1000. Comparison 

of the results for the IRAS- and optically-selected samples shows that there is a population of 
, , galaxies containing a large proportion of cold dust that is unrepresented in the IRAS sample. 

' We derive local submillimetre luminosity and dust mass functions, both directly from our 

optically-selected SLUGS sample, and by extrapolation from the IRAS PSCz survey using 
the method of Serjeant & Harrison (by extrapolating the spectral energy distributions of the 
IRAS PSCz survey galaxies out to 850 |J.m we probe a wider range of luminosities than probed 
directly by the SLUGS samples), and find excellent agreement between the two. We find them 
to be well-fitted by Schechter functions except at the highest luminosities. We find that as a 
consequence of the omission of cold galaxies from the IRAS sample the luminosity function 
presented in our earlier work is too low by a factor of 2, reducing the amount of cosmic 
evolution required between the low-z and high-z Universe. 

Key words: surveys - dust, extinction - galaxies: ISM - galaxies: luminosity function, mass 
function - infrared: galaxies - submillimetre. 



1 INTRODUCTION 

Relatively little is known about the submillimetre properties of 
'normal' galaxies in the local Universe. The advent of IRAS in the 
1980s brought the first investigations of dust in relatively large sam- 
ples of galaxies (e.g. Devereux & Young 1990), yet the limitations 
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of investigating dust at far-IR wavelengths are marked; the strong 
temperature dependence of thermal emission means that even a 
small amount of warm dust can dominate the emission from a sub- 
stantially larger proportion of cold dust, and IRAS is only sensi- 
tive to dust with T > 30 K. IRAS studies of 'normal' galaxies (e.g. 
Devereux & Young 1990) found a high value of the gas-to-dust 
ratio (~1000), an order of magnitude higher than found for the 
Milky Way (~160; Dunne et al. 2000), indicating that IRAS may 
have 'missed' ~90% of the dust in late-type galaxies. IRAS also 
revealed relatively little about the dust in early-type galaxies, since 
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only ~15% of ellipticals were detected by IRAS (Bregman et al. 
1998). 

The next major step in the study of dust in galax- 
ies is to make observations in the submillimetre waveband 
(100 |J.m ^ A ^ 1 mm) since the 90% of dust that is too cold to 
radiate in the far-IR will be producing most of its emission in this 
waveband. The advent of the SCUBA camera on the James Clark 
Maxwell Telescope (JCMT)^ (Holland et al. 1999) opened up the 
submillimetre waveband for astronomy and made it possible, for 
the first time, to investigate the submillimetre emission of a large 
sample of galaxies; prior to SCUBA only a handful of submillime- 
tre measurements had been made of nearby galaxies, using single- 
element bolometers. In particular, in contrast to the extensive sur- 
vey work going on at other wavelengths, prior to SCUBA it was 
not possible to carry out a large survey in the submillimetre wave- 
band. SCUBA has 2 bolometer arrays (850 |i^m and 450 ^.m) which 
operate simultaneously with a field of view of ~2 arcminutes. At 
850 \im SCUBA is sensitive to thermal emission from dust with 
fairly cool temperatures (T > 10 K) so crucially, whereas IRAS 
was only sensitive to warmer dust (T > 30 K), SCUBA should 
trace most of the dust mass. 



1.1 A local submillimetre galaxy survey 

A survey of the dust in nearby galaxies is also important because 
of the need to interpret the results from surveys of the distant Uni- 
verse. Many deep SCUBA surveys have been carried out (Small, 
Ivison & Blain 1997; Hughes et al. 1998; Barger et al. 1998, 1999; 
Blain et al. 1999a; Eales et al. 1999; Lilly et al. 1999; Mortier et 
al. 2005), but studies of the high redshift Universe, and in particu- 
lar studies of cosmological evolution (Eales et al. 1999; Blain et al. 
1999b), have until now depended critically on asstunptions about, 
rather than measurements of, the submillimetre properties of the lo- 
cal Universe. Prior to the existence of a direct local measurement of 
the submillimetre luminosity function (LF) most deep submillime- 
tre investigations have started from a local IRAS 60 iJ.m LF, extrap- 
olating out to submillimetre wavelengths by making assimiptions 
about the average FIR-submm SED. However, as shown by Dunne 
et al. (2000), this underestimates the local submillimetre LF, and 
thus a direct measurement of the local submillimetre LF is vital for 
overcoming this significant limitation in the interpretation of the 
results of high-redshift surveys. 

The ideal method of carrying out a submillimetre survey of 
the local Universe would be to survey a large area of the sky and 
then measure the redshifts of all the submillimetre sources found by 
the survey. However, with current submillimetre instruments such 
a survey is effectively impossible since, for example, the field of 
view of SCUBA is only ~2 arcminutes. The alternative method, 
and the only one that is currently practical, is to carry out targeted 
submillimetre observations of galaxies selected from statistically 
complete samples selected in other wavebands. With an important 
proviso, explained below, it is then possible to produce an unbi- 
ased estimate of the submillimetre LF using 'accessible volume' 
techniques (Avni & Bahcall 1980) (see Section 5). 

To this end several years ago we began the SCUBA Local 
Universe Galaxy Survey (SLUGS). In Papers I and II Duime et 
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al. (2000, hereafter DOO) and Dunne & Eales (2001, hereafter 
DEOl) presented the results of SCUBA observations of a sam- 
ple selected at 60 |im (the //MS-Selected sample, hereafter the IRS 
sample). This paper presents the results of SCUBA observations of 
an optically-selected sample (hereafter the OS sample). 

The accessible volume method will produce unbiased esti- 
mates of the LF provided that no class of galaxy is unrepresented 
in the samples used to construct the LF. DOO produced the first di- 
rect measurements of the submillimetre LF and dust mass func- 
tion (the space-density of galaxies as a function of dust mass) 
using the IRS sample, but this LF would be biased if there ex- 
ists a 'missed' population of submillimetre-emitting galaxies, i.e. 
a population that is not represented at all in the IRS sample. In 
this earlier work we found that the slope of the submillimetre LF 
at lower luminosities was steeper than —2 (a submillimetre 'Gi- 
bers' Paradox'), which indicated that the IRS sample may not be 
fully representative of all submillimetre-emitting sources in the lo- 
cal Universe. This 'missed' population could consist of cold-dust- 
dominated galaxies, i.e. galaxies containing large amounts of 'cold' 
dust (at T < 25 K), which would be strong emitters at 850 \im but 
weak 60 p-m-emitters. The OS sample is selected on the basis of 
the optical emission from the galaxies and, unlike the IRS sample 
which was biased towards warmer dust, the OS sample should be 
free from dust temperature selection effects. The results from the 
OS sample will therefore test the idea that our earUer IRS sample 
LF was an underestimate. 

1.2 Previous investigations of cold dust in galaxies 

The paradigm for dust in galaxies is that there are two main com- 
ponents: (i) a warm component (T > 30 K) arising from dust grains 
near to star-forming regions and heated by young (OB) stars, and 
(ii) a cool 'cirrus' component (T= 15-25 K) arising from diffuse 
dust associated with the HI and heated by the general interstel- 
lar radiation field (ISRF) (Cox, Kriigel & Mezger 1986; Lonsdale 
Persson & Helou 1987; Rowan-Robinson & Crawford 1989). IRAS 
would only have detected the warm component, hence using IRAS 
fluxes alone to estimate dust temperature would result in an over- 
estimate of the dust temperature and an imderestimate of the dust 
mass. Conversely, using the submillimetre to estimate dust masses 
has clear advantages. The flux is more sensitive to the mass of the 
emitting material and less sensitive to temperature in the Rayleigh- 
Jeans part of the Planck function, which is sampled when looking 
at longer submillimetre wavelengths. 

Studies at the longer wavelengths (170-850 \im; e.g. ISO, 
SCUBA) have confirmed the existence of cold dust components 
(15 < Td < 25 K), in line with the theoretical prediction of grain 
heating by the general ISRF (Cox et al. 1986), both in nearby spiral 
galaxies and in more IR-luminous/interacting systems (Guelin et 
al. 1993, 1995; Sievers et al. 1994; Sodroski et al. 1994; Neininger 
et al. 1996; Braine et al. 1997; Dumke et al. 1997; Alton et al. 
1998a,b, 2001; Haas et al. 1998; Davies et al. 1999; Frayer et al. 
1999; Papadopoulos & Seaquist 1999; Xilouris et al. 1999; Haas 
et al. 2000; DEOl; Popescu et al. 2002; Spinoglio et al. 2002; Hip- 
pelein et al. 2003; Stevens, Amure & Gear 2005). Many of these 
authors find an order of magnitude more dust than IRAS observa- 
tions alone would indicate. Alton et al. (1998a), for example, find, 
by comparing their 200 ^m images of nearby galaxies to B-band 
images, that the cold dust has a greater radial extent than the stars, 
and conclude that IRAS 'missed' the majority of dust grains lying 
in the spiral disks. Other studies find evidence of cold dust compo- 
nents in a large proportion of galaxies. Contursi et al. (2001) find 
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evidence of a cold dust component (T ~ 22 K) for most of their 
sample of late-type galaxies; Stickel et al. (2000) find a large frac- 
tion of sources in their 170 iJ,m survey have high 5i7o /Sioo flux 
ratios and suggest this indicates a cold dust component (T < 20 K) 
exists in many galaxies; Popescu et al. (2002) find, for their sample 
of late-type (later than SO) galaxies in the Virgo Cluster, that 30 out 
of 38 galaxies detected in all three observed wavebands (60, 100 
and 170 iJ.m) exhibit a cold dust component. 

An additional property of dust that can be investigated with 
submillimetre measurements is the dust emissivity index 0. Dust 
radiates as a modified Planck function (a 'grey-body'), modified by 
the emissivity term such that Qem <xv^ . Until recently the value of 
P was quite uncertain, with suggested values lying between 1 and 
2 (Hildebrand 1983). Recent multi-wavelength studies of galaxies 
including submillimetre observations, however, have consistently 
fotmd 1.5 ^ /3 ^ 2 with /3=2 tending to be favoured (Chini et al. 
1989; Chini & Krttgel 1993; Braine et al. 1997; Alton et al. 1998b; 
Bianchi et al. 1998; Prayer et al. 1999; DEOl). This agrees with 
the values found in COB£/FIRAS studies of the diffuse ISM in the 
Galaxy (Masi et al. 1995; Reach et al. 1995; Sodroski et al. 1997). 



1.3 The scope of this paper 

This paper presents the results from the SCUBA Local Universe 
Galaxy Survey (SLUGS) optically-selected sample. This OS sam- 
ple is taken from the Center for Astrophysics (Cf A) optical redshift 
survey (Huchra et al. 1983), and includes galaxies drawn from right 
along the Hubble sequence. 

In Section 2 we discuss our observation and data reduction 
techniques. Section 3 presents the sample and the results. Sec- 
tion 4 presents an analysis of the submillimetre properties of the 
sample. In Section 5 we present the local submillimetre lumi- 
nosity and dust mass functions. We assume a Hubble constant 
i7o=75 kms~^ Mpc~^ throughout. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 The sample 

This OS sample is taken from the Center for Astrophysics 
(CfA) optical redshift survey (Huchra et al. 1983), which is 
a magnitude-limited sample of optically-selected galaxies, com- 
plete to ruB ^ 14.5 mag. It has complete information on magni- 
tude, redshift and morphological-type, and also avoids the Galac- 
tic plane. The OS sample consists of all galaxies in the CfA 
sample lying within three arbitrary strips of sky: (i) all decli- 
nations from (B1950.0) 16.1 < RA < 21.5, (ii) all RAs from 
15 < Dec < 16 and (iii) RAs from 9.6 < RA < 12.8 with decU- 
nations from 25 < Dec < 26. We also imposed a lower velocity 
limit of 1900 km s^^ to try to ensure that the galaxies did not have 
an angular diameter larger than the field of view of SCUBA. There 
are 97 galaxies in the CfA survey meeting these selection criteria 
and of these we observed 81 (which were at convenient positions 
given our observing schedule). The OS sample covers an area of 
~ 570 square degrees and is listed in Table 1. Unlike the IRS sam- 
ple which contained many interacting pairs (most of which were 
resolved by SCUBA but not by IRAS), the OS sample contains just 
2 such pairs. 



2.2 Observations 

We observed the OS sample galaxies using the SCUBA bolome- 
ter array at the 15-m James Clark Maxwell Telescope (JCMT) on 
Mauna Kea, Hawaii, between December 1997 and January 2001, 
with a handful of additional observations in February 2003 (due 
to bad data obtained when observed the first time round; see Sec- 
tion 2.3). Observational methods and techniques were similar to 
those for the IRS sample described in DOO. We give a brief descrip- 
tion of these below. 

The SCUBA camera has 2 bolometer arrays (850 ^lm and 
450 nm, with 37 and 91 bolometers respectively) which oper- 
ate simultaneously with a field of view of ~ 2.3 arcminutes at 
850 |im (sUghtly smaller at 450 |im). Beamsizes are measured to 
be ~15 arcsec at 850 ^m and ~8 arcsec at 450 |a.m. Our observa- 
tions were made in 'jiggle-map' mode which, for sources smaller 
than the field of view, is the most efficient mapping mode. Since 
the arrangement of the bolometers is such that the sky is instanta- 
neously undersampled, and since we observed using both arrays, 
the secondary mirror was stepped in a 64-point jiggle pattern in or- 
der to fully sample the sky. The cancellation of rapid sky variations 
is provided by the telescope's chopping secondary mirror, operating 
at 7.8 Hz. Linear sky gradients and the gradual increase or decrease 
in sky brightness are compensated for by nodding the telescope 
to the 'off position every 16 seconds. We used a chop throw of 
120 arcsec in azimuth, except where the galaxy had a nearby com- 
panion, in which case we used a chop direction which avoided the 
companion. 

The zenith opacity r was measured by performing regular 
skydips. The observations were carried out under a wide range 
of weather conditions, with opacities at 850 \im rgso ranging from 
0.12 to 0.52. This means that some galaxies were observed in excel- 
lent conditions (rsso < 0.2) while others were observed in far less 
than ideal conditions. As a result we obtained useful 450 iim data 
for only a fraction of our galaxies. This is discussed in more de- 
tail in Section 2.5. Our observations were centred on the coordi- 
nates taken from the NASA/IPAC Extragalactic Database (NED). 
We made regular checks on the pointing and found it to be gener- 
ally good to ~2 arcsec. The integration times depended on source 
strength and weather conditions. Since most of the OS sample are 
relatively faint submillimetre sources we typically used ~12 in- 
tegrations (~30 mins), although many sources were observed in 
poorer weather and so required longer integration times. 

We calibrated our data by making jiggle maps of Uranus and 
Mars, or, when these planets were unavailable, of the secondary 
calibrators CRL 618 and HL Tau. We took the planet fluxes from 
the JCMT FLUXES program, and CRL 618 and HL Tau were as- 
sumed to have fluxes of 4.56 and 2.32 Jy beam~^ respectively at 
850 [im. 



2.3 Data reduction 

The 850 n-m and 450 ^m data was reduced using the standard 
SCUBA specific tasks in the SURF package (Jenness & Lightfoot 
1998, 2000; Jenness et al. 2002), where possible via the XORACDR 
automated data reduction pipeline (Economou et al. 2004). The 
off-nod position was subtracted from the on-nod in the raw beam- 
switched data and the data was then flat-fielded and corrected for 
atmospheric extinction. 

In order to correct SCUBA data for atmospheric extinction we 
must accurately know the value of the zenith sky opacity, r. Al- 
though less crucial at 850 [im if the observation is made in good 
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Table 1. 850 (jm flux densities and isothermal SED parameters. (Notes on individual objects are listed in Section 3.1). 
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4-11 57 51 

T J 1 J / J 1 


049 n 

7tZD 


n fi4 

U.Dt 


1 fi7 
l.D / 


n n59 

D.DjZ 


019 

U.UIZ 


11 4 
J j.t 


8 

D.O 


J 


NPP 6104 


1 fi 1 10 fi 

lU lU Jv.U 


+35 42 29 


8428 


0.50 


1.76 


<0.033 










IC 1211 


Ifi 1^^ SI 

ID 10 J i-.y 


4-51 no 99 

TJJ vU ZZ 


5fil8 

jDiO 


19 

<s..D. iZ 


^0 51 

<s..D. J J 


n n98 

D.DZO 


000 






-J 


TTPP 1 0^9 '^^ 8 


1 (i 1 7 10 (i 


4-46 05 10 
t'tu vj JV 


5691 


1.57 


3.72 


0.041 


0.009 


31.0 


1.4 


lOM 


IN OL. D 1 z / 


1 A 1 Q 1 1 ^ 
10 IV 1 1 . J 


+j / Jy Uj 


to J i 


^0 1 n 




U.UoD 


u.uzu 






c 
-J 


NGC6120 


16 19 48.1 


+37 46 28 


9170 


3.99 


8.03 


0.065 


0.011 


32.2 


1.5 


8 


NGC6126 


16 21 27.9 


+36 22 36 


9759 


<0.15 


<0.43 


0.023 


0.008 






-2 


NGC6131 


16 21 52.2 


+38 55 57 


5117 


0.72 


2.42 


0.054 


0.013 


28.6 


1.2 


6 


NGC 6137 


16 23 03.1 


+37 55 21 


9303 


<0.18 


<0.53 


0.029 


0.010 






-5 


NGC 6146 


16 25 10.3 


+40 53 34 


8820 


<0.12 


<0.48 


0.028 


0.007 






-5 


NGC 6154 


16 25 30.4 


+49 50 25 


6015 


<0.15 


<0.36 


<0.040 








1 


NGC 6155 


16 26 08.3 


+48 22 01 


2418 


1.90 


5.45 


0.116 


0.022 


29.8 


1.2 


6 


UGC 10407 


16 28 28.1 


+41 13 05 


8446 


1.62 


3.12 


0.026 


0.009 


32.8 


1.5 


lOM 


NGC 6166 


16 28 38.4 


+39 33 06 


9100 


^0.10 


''0.63 


0.073 


0.017 


26.2 


0.6 


-5 


NGC 6173 


16 29 44.8 


+40 48 42 


8784 


<0.17 


<0.23 


<0.024 








-5 


NGC 6189 


16 3140.9 


+59 37 34 


5638 


0.75 


2.57 


0.072 


0.019 


28.6 


1.1 


6 


NGC 6190 


16 32 06.7 


+58 26 20 


3351 


0.58 


2.37 


0.099 


0.024 


28.0 


1.0 


6 



The SCUBA Local Universe Galaxy Survey - 



Table 1 - continued 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


Name 


R.A. 


Decl. 


cz 




Sioo 




csso 




P 


Type 




(J2000) 


(J2000) 


/I — 1\ 
(kms ) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(K) 






NGC 6185 


16 33 17.8 


+35 20 32 


10301 


0.17 


0.56 


<0.030 








1 


UGC 10486 


16 37 34.3 


+50 20 44 


6085 


<0.19 


<0.60 


<0.029 








-3 


NGC 6196 


16 37 53.9 


+36 04 23 


9424 


<0.12 


<0.44 


<0.023 








-3 


UGC 1050{) 


16 38 59.3 


+57 43 27 


5218 


**0.16 


■^*0.71 


<0.028 











IC 5090 


21 11 30.4 


-02 01 57 


9340 


3.04 


7.39 


0.118 


0.017 


31.6 


1.2 


1 


IC 1368 


21 14 12.5 


+02 10 41 


3912 


4.03 


5.80 


0.047 


0.011 


37.6 


1.3 


1 


NGC 7047 


21 16 27.6 


-00 49 35 


5626 


0.43 


1.65 


0.055 


0.013 


28.0 


1.1 


3 


NGC 7081 


21 31 24.1 


+02 29 29 


3273 


1.79 


3.87 


0.044 


0.010 


32.8 


1.3 


3 


NGC 7280 


22 26 27.5 


+16 08 54 


1844 


<0.12 


<0.48 


<0.040 








-1 


NGC 7442 


22 59 26.5 


+15 32 54 


7268 


0.78 


2.22 


0.046 


0.009 


31.0 


1.1 


5 


NGC 7448t 


23 00 03.6 


+15 58 49 


2194 


7.23 


17.43 


0.193 


0.032 


31.0 


1.4 


5 


NGC 7461 


23 01 48.3 


+ 15 34 57 


4272 


<0.176 


<0.64 


<0.022 








-2 


NGC 7463 


23 01 51.9 


+15 58 55 


2341 


U 


U 


0.045 


0.010 






3M 


111 ZW 093 


23 07 21.0 


+15 51 11 


14962 


0.48 


<3.16 


<0.026 








lOZ 


111 ZW 095 


23 12 43.3 


+15 54 12 


7506 


<0.09 


<0.80 


<0.019 








lOZ 


UGC 12519 


23 20 02.7 


+15 57 10 


4378 


0.76 


2.59 


0.074 


0.016 


29.2 


1.1 


5 


NGC 7653 


23 24 49.3 


+15 16 32 


4265 


1.31 


4.46 


0.112 


0.020 


28.6 


1.2 


3 


NGC 7691 


23 32 24.4 


+15 50 52 


4041 


0.53 


1.67 


<0.025 








4 


NGC 7711 


23 35 39.3 


+15 18 07 


4057 


<0.15 


<0.50 


<0.027 








-2 


NGC 7722 


23 38 41.2 


+15 57 17 


4026 


0.78 


3.03 


0.061 


0.015 


26.8 


1.4 






(1) Most commonly used name. 

(2) Right ascension, J2000 epoch. 

(3) Declination, J2000 epoch. 

(4) Recessional velocity taken from NED. [The NASA/IPAC Extragalactic Database (NED) is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration.] 

(5) 60 \im flux from the IRAS Faint Source Catalogue (Moshir et al. 1990); upper limits listed are measured using SCANPI as 

described in Section 2.8. 

(6) 100 i^Ti flux from the IRAS Faint Source Catalogue (Moshir et al. 1990); upper limits Usted are measured using SCANPI as 
described in Section 2.8. 

(7) 850 Mm flux (this work). 

(8) Error on 850 |jm flux, calculated as described in Section 2.6. 

(9) Dust temperature derived from a single-component fit to the 60, 100 and 850 |im data points, as described in Section 3.2. 

(10) Emissivity index derived from the single-component fit, as described in Section 3.2. 

(11) Hubble type (t-type) taken from the LEDA database; we have assigned t=10 to any multiple systems umesolved by IRAS 
or SCUBA (indicated by 'lOM') and any systems with no type listed in LEDA (indicated by 'lOZ'; these 2 objects are listed as 
'compact' sources in NED); all other types marked 'M' are hsted as multiple systems in LEDA. 

^ Part of a close or interacting pair which was resolved by SCUBA. Fluxes here are the individual galaxy fluxes; fluxes measured 
for the combined pair are given in Table 2. 
U Unresolved by IRAS. 

* The IRAS flux is our own SCANPI measurement (see Section 2.8); any individual comments are listed in Section 3.1. 

* SCANPI measurements and fitted values should be used with caution (see Section 2.8). 
§ The coordinates of this object refer to one galaxy (NEDOl) of a the pair UGC 10325. 

t Objects are also in the Paper I ffiAS-selected sample (DEOO). 



Table 2. Combined SCUBA fluxes for pairs umesolved by IRAS. (Notes on individual objects are hsted in Section 3.1). 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


Name 


R.A. 


Decl. 


cz 


560 


Sioo 


SsBO 


CSBO 


rp 

dust 


P 


Type 




(J2000) 


(J2000) 


(kms-i) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(K) 






NGC 3799/3800 


11 40 11.4 


+15 20 05 


3312 


4.81 


11.85 


0.135 


0.035 


29.8 


1.5 


lOM 


NGC 5953/4 


15 34 33.7 


+15 11 49 


1966 


10.04 


18.97 


0.273 


0.034 


35.2 


1.1 


lOM 



Note. Columns have the same meanings as in Table 1. 
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weather (tsso <0.3) and at low airmass, in worse weather or at 
450 |im the measured source flux can be severely affected by an er- 
ror in T. r is most commonly estimated either by performing a sky- 
dip or by extrapolating to the required wavelength (using relations 
given in the JCMT literature and in Archibald et al. (2002)) from 
polynomial fits to the continuous measurements of t at 225GHz 
made at the nearby Caltech Submillimetre Observatory. Since sky- 
dips are measured relatively infrequently, the polynomial fits to the 
CSO T225 data are recommended in the JCMT literature to be the 
more reliable way of estimating r for both SCUBA arrays. As 
such, for both 850 |i^m and 450 ^.m data we have wherever pos- 
sible (the large majority of observations) used the derived CSO 
opacity at 225GHz {tcso)- Where Tcso values were not available 
the opacities were derived from 850 \im skydip measurements (at 
450 (i.m using the tsso -to-T45() relation described in the JCMT lit- 
erature and Archibald et al. (2002)). 

Noisy bolometers were noted but not removed at this stage (it 
was frequently found to be the case that flagging a noisy bolometer 
as 'bad' creates even worse noise spikes in the final map around 
the position of the removed bolometer data). Large spikes were re- 
moved from the data using standard SURF programs. 

The nodding and chopping should remove any noise which 
is correlated between the different bolometers. In reality, since the 
data was not observed in the driest and most stable conditions the 
signal on different bolometers was often highly correlated due to 
incomplete sky subtraction. In the majority of cases we used the 
SURF task REMSKY, which takes a set of user-specified bolometers 
to estimate the sky variation as a function of time. More explic- 
itly, in each time step REMSKY takes the median signal from the 
specified sky bolometers and subtracts it from the whole array. To 
ensure that the sky bolometers specified were looking at sky alone 
and did not contain any source emission we used a rough SCUBA 
map together with optical {Digitised Sky Survey'^ (DSS)) images as 
a guide when choosing the bolometers, though in this sample there 
are so few bright sources that in the majority of cases all bolometers 
could be safely used. 

Even after this step, however, due to the relatively poor con- 
ditions in which much of the data was observed the residual sky 
level was sometimes found to vary Unearly across the array, giv- 
ing a 'tilted plane' on the array. Moreover, in a number of cases a 
noisy 'striped' sky (due possibly to some short-term instrumenta- 
tion problem) was found. Though the SURF task REMSKY was de- 
signed to remove the sky noise, it is relatively simplistic and cannot 
remove such spatially varying 'tilted' or 'striped' sky backgrounds. 
In these cases, as for the IRS sample (DOO), we used one of our own 
programs in place of REMSKY. In a handful of cases the 'striped' 
sky was so severe that it could not be removed, so these objects 
were re-observed in February 2003. 

Once the effects of the sky were removed the data was de- 
spiked again and the final map produced by re-gridding the data 
into a pixel grid to form an image on 1 arcsecond pixels. Where 
there were multiple data-sets for a given source they were binned 
together into a co-added final map. In these cases each data set 
was weighted prior to co-adding using the SURF task SETBOLWT, 
which calculates the standard deviation for each bolometer and then 

The Digitised Sky Surveys were produced at the Space Telescope Science 
Institute under U.S. Government grant NAGW-2166. The images of these 
surveys are based on photographic data obtained using the Oschin Schmidt 
Telescope on Palomar Mountain and the UK Schmidt Telescope. The plates 
were processed into the present digital form with the permission of these 
institutions. 



calculates weights relative to the reference bolometer (the central 
bolometer in the first input map). This method is therefore only 
suitable if there are no very bright sources present in the central 
bolometer (if a bright source was present we weighted each dataset 
using the inverse square of its measured average noise). This step 
also ensures that noisy bolometers contribute to the final map with 
their correct statistical weight. 

2.4 850 ^m flux measurement 

The fluxes were measured from the SCUBA maps by choosing 
a source aperture over which to integrate the flux, such that the 
signal-to-noise was maximised. The extent of the galaxy in the op- 
tical (DSS) images and the extent of the submillimetre source on 
the S/N map (see Section 2.7) were used to select an aperture that 
included as much of the submillimetre flux of the galaxy as pos- 
sible while minimising the amount of sky included. Note, the op- 
tical images in Figure 1 are shown stretched for optimum contrast 
- however, apertures for flux measurement were drawn for a more 
modest optical extent, as seen at a standard level of contrast. 

Conversion of the measured aperture flux in volts to Janskys 
was carried out by measuring the calibrator flux for that night using 
the same aperture as for the object. The orientation of the aperture 
(relative to the chop throw) was also kept the same as for the object, 
as particularly for more elliptical apertures this has a significant 
effect. 

Objects are said to be detected at > 3c7 if either: (a) the peak 
S/N in the S/N map was > 3cr or (b) the flux in the aperture was 
greater than 3 times the noise in that aperture (where the noise is 
defined as described in Section 2.6). 

2.5 450 ^m data 

Due to the increased sensitivity to weather conditions at 450 |im, 

sources emitting at 450 |J.m will only be detected if they are rel- 
atively bright at 850 ]xm. This, together with the wide range of 
observing conditions for this sample, meant that we found useful 
450 ]xm data for only 19 objects. 

Where possible the 450 ^.m emission was measured in an 
aperture the same size as used for the 850 |im data. In some cases 
a smaller aperture had to be used for the 450 [im data, and these 
individual cases are discussed in Section 3.1. 



2.6 Error analysis 

The error on the flux measurement is made up of three components: 

• A background sky subtraction error cTsky due to the uncer- 
tainty in the sky level. 

• A shot (Poisson) noise term ashot due to pixel-to-pixel varia- 
tions within the sky aperture. Unlike CCD images, in SCUBA maps 
the signal in adjoining pixels is correlated; this correlated noise de- 
pends on a number of factors, including the method by which the 
data is binned at the data reduction stage. This has been discussed 
in some detail by DOO, who find that a correction factor is required 
for each array to account for the fact that pixels are correlated; they 
find the factor to be 8 at 850 ^.m and 4.4 at 450 ^.m. 

• A calibration error term acai which for SCUBA observations 
at 850 |j,m is typically less than 10%. We have therefore assumed a 
conservative caUbration factor of 10% at 850 ia.m. The calibration 
error at 450 ia.m was taken to be 15%, following DEOl. 
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Table 3. 450 \na flux densities and two-component SED parameters. (Notes on individual objects 
are listed in Section 3.1). 



(1) 


(2) 






P) 


(o) 


(7) 


(8) 


(9) 


Name 




O"450 


3 

■5850 




Tc 






Lfir 




(Jy) 


(Jy) 




(K) 


(K) 


(log Mo) 


(log Lq) 




0.944 


0.236 


17.18 


34 


18 


37 




10.33 


NGC 99 


0.490 


0.182 


7.73 


47 


17 


542 


1.12 


10.08 


Nfir 803 1 


0.631 


0.196 


6.79 


33 


18 


92 


1.02 


9.46 


NGC 3689* 


1.045 


0.357 


10.30 


59 


23 


910 


7.13 


10.16 


PGC 3'i9'i2 


0.421 


0.116 


8.26 


58 


18 


1859 


7.31 


9.77 


NGC 3987 


1.110 


0.319 


5.98 


44 


22 


279 


7.85 


10.78 


IC 979§ 


0.874 


0.341 


15.39 












NGC 5953/4 


2.879 


0.683 


10.54 


54 


21 


277 


7.33 


10.28 


NGC 5980 


1.398 


0.495 


5.53 


43 


18 


321 


8.06 


10.53 


NGC 6090 


0.803 


0.180 


8.82 


55 


22 


122 


8.09 


11.29 


NGC 6120 


0.528 


0.127 


8.08 


45 


24 


76 


7.96 


11.17 


NGC 6155* 


0.381 


0.135 


3.30 


30 


20 


7 


6.92 


9.80 


NGC 6190* 


0.880 


0.308 


8.89 


56 


18 


2684 


7.16 


9.85 


IC 5090 


1.018 


0.240 


8.66 


52 


21 


346 


8.28 


11.19 


IC 1368 


0.425 


0.137 


9.10 


55 


23 


110 


7.10 


10.37 


NGC 7081 


0.241 


0.067 


5.43 


32 


20 


6 


6.98 


9.93 


NGC 7442 


0.410 


0.099 


9.02 


54 


20 


665 


7.70 


10.45 


UGC 12519 


0.408 


0.108 


5.54 


28 


17 


12 


7.57 


10.02 


NGC 7722 


0.595 


0.148 


9.78 


54 


20 


1224 


7.33 


10.04 



(1) Most commonly used name. 

(2) 450 yum flux (this work). 

(3) Error on 450 [jm flux, calculated as described in Section 2.6. 

(4) Ratio of 450- to 850- fluxes. 

(5) Warm temperature using /? = 2. 

(6) Cold temperature using /3 = 2. 

(7) Ratio of cold-to-warm dust. 

(8) Dust mass calculated using parameters in columns (5)-(7). 

(9) FIR luminosity (40-1000(jm) integrated under the two-component SED. 
* Some caution is advised (see Section 3.1). 

§ The data could not be fitted with a 2-component model (see Section 3.1). 

t Not well-fitted by two-component model using 850 \xm data point; fitted parameters here are 
from 2-component fit to the 60, 100, 450 and 170 |am (ISO) data points; (see Section 3.1). 
tThe data for NGC 803 are also weU-fitted by the parameters: ru,=60K, Tc=19K, -^=2597, 
log Md2=6.99 and log L / j^=9.48, (see Section 3.1). 



The relationships used to calculate the noise terms are as fol- 
lows: 

and 

CTshot ~ Sapixy^ Nap or (Tshot = 4:-'io'pix\/ Nap 

for 850 |a.m and 450 \un flux measurements respectively. The error 
in the mean sky a-ms = S.D./s/n, where S.D. is the standard de- 
viation of the mean sky values in n apertures placed on off-source 
regions of the map. Nap is the number of pixels in the object aper- 
ture; (jpix is the mean standard deviation of the pixels within the 
sky apertures. The total error for each flux measurement is then 
given by 

/ 2 . 2 I 2 -.1/2 /1^ 
Ctot = {CTaky + ""shot + '^cal) W 

as for the IRS sample. This error analysis is discussed in detail in 
DOOand DEOl. 

850 ^.m fluxes were found to have total errors atot typically in 
the range 15-30 %. 450 jxm fluxes were found to have total errors 
ctot typically in the range 25-35 %. Note, the atot used to deter- 



mine whether a source was detected at the 3(t level is defined as in 
Equation 1 but without the calibration error term. 



2.7 S/N maps 

Unlike the IRS sources the OS sources were not selected on the 
basis of their dust content. Many of the OS sources, especially the 
early types, are close to the limit of detection. Also, it is often hard 
to assess whether a source is detected, or whether some feature of 
the source is real, due to the variability of the noise across the ar- 
ray. This is due both to an increase in the noise towards the edge 
of each map, caused by a decrease in the number of bolometers 
sampling each sky point, and to individual noisy bolometers. For 
this reason we used the method described in DOO to generate ar- 
tificial noisemaps, wiiich we used with our real maps to produce 
signal-to-noise maps. The real maps and the artificial maps were 
first smoothed (using a 12 pixel FWHM) before creating the S/N 
map. 

We used these S/N maps to aid in choosing the aperture for 
measuring the 850 [xm flux (Section 2.4). We have also presented 
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S/N maps of each source (see Section 3), as this makes it easier to 
assess the reality of any features in the maps. 

2.8 IRAS fluxes 

IRAS 100 (im and 60 (i.m fluxes, where available, were taken from 
the IRAS Faint Source Catalogue (Moshir et al. 1990; hereafter 
FSC) via the NED database. Where literature fluxes were unavail- 
able die NASA/IPAC Infrared Science Archive (IRSA) SCANPl 
(previously ADDSCAN) scan coadd tool was used to measure a 
flux from the IRAS survey data. 

The small number of SCANPl fluxes are indicated by 's' in 
Table 1, and any special cases are discussed individually in Sec- 
tion 3.1. We take SCANPl fluxes to be detections if the measure- 
ments are formal detections at > 4.5cr at 100 i^m or > 4cr at 60 [im, 
which Cox et al. (1995) conclude are actually detections at the 98% 
confidence level. Otherwise we give a 98% confidence upper limit 
(4.5a at 100 |j.m or 4a at 60 ]im) using the la error found from 
SCANPl (again following Cox et al. (1995)). If both fluxes are 
SCANPl measurements we mark the subsequent fltted values by 
'*' if there is any doubt as to their viability (for example possible 
source confusion, confusion with galactic cirrus, or no literature 
IRAS fluxes in NED for either band). 



3 RESULTS 

We detected 52 of the 81 galaxies in the OS sample. Table 1 lists 
the 850 ^m fluxes and other parameters. For interacting systems 
resolved by SCUBA but not resolved by IRAS the 850 \im fluxes 
given are for the individual galaxies; the 850 |j.m fluxes measured 
for the combined system are listed in Table 2 along with the IRAS 
fluxes. Table 3 lists the 450 (xm fluxes for the 19 galaxies which 
are also detected at the shorter wavelength. The galaxies detected 
in the OS sample are shown in Figure 1, with our 850 |im SCUBA 
S/N maps overlaid onto optical (DSS) images. Comments on the 
individual maps are given in Section 3.1. 

The 850 |im images have several common features. Firstly, 
we find that many spiral galaxies exhibit two peaks of 
850 ^.m emission, seemingly coincident with the spiral arms. This 
is most obvious for the more face-on galaxies (for example NGC 99 
and NGC 7442), but it is also seen for more edge-on spirals (e.g. 
NGC 7047 and UGC 12519). This 'two-peak' morphology is not 
seen for all the spirals, however. Some, for example NGC 3689, are 
core-dominated and exhibit a single central peak of submillimetre 
emission, while others (NGC 6131 and NGC 6189 are clear ex- 
amples) exhibit a combination of these features, with both a bright 
nucleus and peaks coincident with the spiral arms. In a number of 
cases the 850 |J.m peaks clearly follow a prominent dust lane (e.g. 
NGC 3987 and NGC 7722). These results are consistent with the 
results of numerous mm/submm studies. For example, Sievers et 
al. (1994) observe 3 distinct peaks in NGC 3627 and note that the 
two outer peaks are coincident with the transition region between 
the central bulge and the spiral arms - they also observe dust emis- 
sion tracing the dust lanes of the spiral arms; Guelin et al. (1995), 
Bianchi et al. (2000), Hippelein et al. (2003) and Meijerink et al. 
(2005) observe a bright nucleus together with extended dust emis- 
sion tracing the spiral arms. Many of the features seen in our OS 
sample 850 ]im maps are also found by Stevens et al. (2005) in their 
SCUBA observations of nearby spirals. 

Secondly, we find that a number of galaxies appear to be ex- 
tended at 850 ia.m compared to the optical emission seen in the DSS 



images. In many cases this extended 850 iJ.m emission appears to 
correspond to very faint optical features, as can be seen for NGC 
7081 and NGC 7442 in Figure 1. In order to investigate this fur- 
ther we have already carried out follow-up optical imaging for ~ 
half the sample detected at 850 |j.m, to obtain deeper images than 
available from the DSS. The results and discussion of this deeper 
optical data will be the subject of a separate paper (Vlahakis et al., 
in preparation). 



3.1 Notes on individual objects 

In the following discussion of individual objects we note that since 
the number of bolometers sampling each sky point decreases to- 
wards the edges of the submillimetre maps the noise increases 
towards the edge of the maps. Although the S/N maps in Fig- 
ure 1 were produced using artificial noisemaps (Section 2.7) which 
should normally account for this effect there are certain circum- 
stances, such as a 'tilted-sky' (see Section 2.3) or the very noisiest 
bolometers, where residual 'noisy' features may remain in the S/N 
maps. This means that any submillimetre emission in Figure 1 seen 
beyond the main optical extent (and away from the centre of the 
map) should be regarded with some caution. However, in order to 
aid distinction between probable residual noisy features in the S/N 
map and potential extended submillimetre emission we have made 
a thorough investigation of each individual map. In the following 
discussion, unless otherwise stated we have found all ^ 2cr submil- 
limetre peaks away from the main optical galaxy to be associated 
with noisy bolometers or a tilted sky. 

UGC 148. Data points for this object are not well-fitted by 
the two-component dust model (Section 3.2), probably due to the 
850 i^m flux having been underestimated - the 850 ^.m S/N con- 
tours shown in Figure 1 for this object show evidence of a residual 
tilted sky plane (the sky is more positive on one side of the map 
than the other), suggesting that sky removal techniques may have 
been inadequate in this case and that therefore the source flux may 
have been under- (or over-) estimated. Also, since the E-NE part of 
the galaxy is coincident with noisy bolometers in the 850 |j.m map 
the flux-measurement aperture was drawn to avoid this region, so 
the 850 i^m flux may be vmderestimated. However, an additional 
data point at 170 |J.m (ISO) is available from the literature (Stickel 
et al. 2000, 2004). Using the 60, 100, 170 and 450 ^m fluxes we 
find that the data points are well-fitted by the two-component dust 
model (we take an average of all 170 \xm fluxes available, see Sec- 
tion 3.2), and these results are listed in Table 3 and the SED is 
shown in Figure 2. 

NGC 99. The submillimetre emission follows the spiral arm 
structure. The 2a peak to the NE of the galaxy is not associated 
with any noisy bolometers. 

NGC 786. The submillimetre emission to the SE of the galaxy 
is not associated with any noisy bolometers but we note that this 
object was observed in very poor weather. 

NGC 803. None of the submiUimetre peaks are associated 
with noisy bolometers. Due to the high ratio of S25/S60 good 
two-component SED fits (Section 3.2) to the 4 data points (60, 100, 
450 and 850 |j.m) can be achieved for two quite different values of 
the warm component temperature. In addition to the parameters 
listed in Table 3 a good fit is also found with the foUowing 
parameters: T„=60K, Tc=19K, ^=2597, log Md2=6.99 and 
logL/„=9.48. 

UGC 5342. This observation had a tilted sky. 

NGC 3270. All the submillimetre emission away from the 
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main optical extent, and the emission to the N of the galaxy, is 
associated with noisy bolometers. This observation also suffered 
from a tilted sky, which may explain much of the submillimetre 
emission in the N part of the map. However, the emission towards 
the centre of the map, coinciding with the main optical galaxy, is 
not associated with any noisy bolometers. This emission seems to 
occur where the galaxy bulge ends and the inter-arm region begins, 
as was found for NGC 3627 by Sievers et al. (1994). 

NGC 3689. Very few scans were available via SCANPl, so 
IRAS fluxes (and corresponding fitted parameters) for this object 
should be used with caution. 

PGC 35952. The submillimetre emission to the S and SW is 
not associated with any noisy bolometers. The submillimetre peaks 
coincident with the main optical extent of the galaxy appear to fol- 
low the spiral arm structure, as for NGC 6131, and it is possible 
the extended submillimetre emission to the S relates to the very 
extended faint spiral arms seen in the optical. 

NGC 3799/3800. NGC 3799 and NGC 3800 were observed 
in separate maps. NGC 3799 individually is not detected at the 3a 
level (although we measure flux at the 2a level). For NGC 3800 
(shown in Figure 1) most of the integrations for the bolometers to 
the S of the map were unusable, and consequently this part of the 
map is very much noisier. Generally this observation is very noisy, 
and especially bad at 450 ^.m; thus no 450 (im flux is available. In 
fact only the main region of submillimetre emission at the centre 
of the map is in an area free from noisy bolometers, and it is this 
region over which we have measured the submillimetre flux of the 
galaxy. 

The Ssso Usted for NGC 3799/3800 in Table 2 is a conserva- 
tive measurement of the 850 iJ,m emission from the system, the sum 
of the separately measured fluxes from each of the two component 
galaxies. In coadding the two maps there appears to be a 'bridge' of 
850 ]xm emission between the two galaxies, consistent with emis- 
sion seen in the optical (NGC 3799 is to the SW of NGC 3800 in 
Figure 1). However, since this region of the map has several noisy 
bolometers we only measure fluxes for the main optical extent of 
the galaxies. 

NGC 3815. The submiUimetre emission to the NE and W of 



the galaxy is associated with noisy bolometers. However, the arm- 
like submillimetre structures seen extending from the galaxy to the 
N and S arc not. Both of these 'arms' extend in the direction of faint 
optical features seen in the DSS and 2MASS images. The optical 
images also show evidence of extended optical emission around 
the main optical extent (just visible to the NE in Figure 1), but 
2MASS (JHK) images show a band of emission stretching E-W 
between two nearby galaxies on either side of NGC 3815. It seems 
clear that some kind of interaction is taking place in this system, 
and therefore it is perhaps not unlikely that there might also be 
significantly extended submillimetre emission. 

NGC 3920. The submillimetre emission to the E and S is not 
associated with any noisy bolometers. 

NGC 3987. This edge-on galaxy has a prominent dust lane 
in the optical. Though the submillimetre emission follows the dust 
lane it is seen slightly offset. A similar result was found for another 
edge-on spiral by Stevens et al. (2005), who conclude this effect is 
simply an effect of the inclination of the galaxy on the sky. 

NGC 3997. The FSC gives an upper hmit <3.101 Jy for Sioo, 
likely due to possible source confusion with NGC 3993. The Sioo 
we measure with SCANPl should therefore be used with some cau- 
tion. 

NGC 4005. This object is detected at 850 [im at only the 2.5a 
level. It is unresolved by IRAS and may be confused with IRAS 
source IRASFl 1554+2524 (NGC 4000). 

UGC 7115. With the exception of the submillimetre emission 
to the SE, none of the submillimetre emission in this map is as- 
sociated with noisy bolometers. However, we found this SCUBA 
observation to have a tilted sky. We estimate that as much as 80% 
of the 850 (im flux from this elliptical may be due to synchrotron 
contamination from a radio source associated with the galaxy (Vla- 
hakis et al., in prep.). 

IC 800. The submillimetre emission to the E of this galaxy 
corresponds to a region of the map which is only slightly noisy, 
making it unlikely any residual features of this noise remains in the 
S/N map (Section 2.7). However, we also note that this object was 
observed in poor weather. 

NGC 4712. None of the submilUmetre emission is associated 
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with noisy bolometers, witti tiie exception of the 2a peak closest to 
the galaxy in the arm-like structure extending to the E. However, 
we note that this arm-like feature, though faint, is also seen in the 
optical (though not reproduced in the optical image in Figure 1); 
it appears to originate from the main galaxy extent, where there 
appears to be a significant amount of dust obscuration. 

UGC 8902. The region of submillimetre emission to the E, far 
SW and far S are all associated with noisy bolometers. However, 
the 4a submillimetre peak lying to the S/SE beyond the main opti- 
cal extent, at a similar declination to the small galaxy to the SE, is 
not associated with any noisy bolometers. This submilUmetre emis- 
sion is consistent with the fact that the overall emission associated 
with the galaxy is offset to the S/SE with respect to the optical. We 
also note that this region in the optical contains a number of faint 
condensations in the direction of the small galaxy to the SE. 

IC 979. Although this galaxy is detected with relatively low 
S/N at 850 (xm it is also detected at 450 \im. We allocate a higher 
450 iJ.m calibration error (25%) for this source, since there were 
no good 450 ^.m calibrator observations that night (calibration was 
achieved by taking the mean results from a number of calibra- 
tors observed that and the previous night). Note also that no two- 
component fit could be made to the data since the 450 \xm data 
point is higher than the 100 iJ.m value, possibly due to the problems 
with calibrating the 450 \i.m data but more likely due to an under- 
estimate of the 100 li^m flux. We note this object was observed in 
poor weather. 

UGC 9110. There appears to be flux present at both 
850 (i.m and 450 (im at the 2cr level, but the maps are very noisy and 
several integrations unusable, most likely due to unstable and dete- 
riorating weather conditions during the observation. This object is 
unresolved by IRAS: IRAS source IRASF141 19-1-1551 (FSC fluxes 
Sioo=2.341 and S6o=0.802 Jy) is likely a blend of UGC9110 and 
companion CGCG103-124 (Condon, Cotton & Broderick 2002). 

NGC 5522. The 850 [im emission to the SE of NGC 5522 is 
associated with a region of the map which is slightly noisy and 
where there are a number of spikes in the data. We note that this 
observation was carried out in poor weather. 

NGC 5953/4 is also in the IRS SLUGS sample. While DOO 
used colour-corrected IRAS fluxes as listed in the IRAS BGS (from 
Soifer et al. (1989)) we present here, as for all the OS sample, fluxes 
from the IRAS FSC. 

NGC 5980. This observation suffered from a tilted sky, which 
potentially explains the extended 850 |j.m emission to the W of the 
galaxy. The aperture used to measure the 450 iJ-m flux was smaller 
than at 850 [im in order to avoid a noisy bolometer, thus source flux 
at 450 |j.m may be tmderestimated. 

IC 1174. This source is just detected at 850 \im at the 3a level. 

UGC 10205. The 850 |J.m map in Figure 1 is a coadd of two 
observations. Since emission at the optical galaxy position is very 
clear in one observation (both at 850 |xm and 450 |J.m) and not in 
the other observation, and since we find no explanation for this, we 
simply coadd the two observations (Section 2.2). The submillimetre 
emission coincident with the main optical galaxy extent, and also 
the peak to the S, are not associated with any noisy bolometers. 
Peaks to the N and W of the galaxy lie in a region of the map which 
is slightly noisy. We note that this observation was carried out in 
less than ideal weather. 

NGC 6090 is a closely interacting/merging pair, and is also in 
the IRS SLUGS sample. 

NGC 6103. The submillimetre emission to the S of the galaxy 
is not associated with any noisy bolometers. This region contains a 



number of faint features seen in optical (DSS) and 2MASS images. 
We note that this object was observed in less than ideal weather. 

IC 1211. We find in the literature no known radio sources as- 
sociated with this elliptical galaxy (NVSS 1 .4GHz 3a upper limit is 
< 1.2 mJy), and therefore cannot attribute the 850 \im flux detected 
here to contamination from synchrotron radiation. 

UGC 10325 (NEDOl) is one galaxy of the pair UGC 10325. 
The SCUBA map is centred on this galaxy (NEDOl), but NED02 
can be seen at the SE edge of the DSS image in Figure 1. Thus all 
fluxes given are for the individual galaxy UGC 10325 NEDOl. 

NGC 6127. We find in the literature no known radio sources 
associated with this elliptical galaxy (NVSS 1.4GHz 3a upper limit 
is <1.2mJy), and therefore cannot attribute the 850 ji-m flux de- 
tected here to contamination from synchrotron radiation. The 4a 
submillimetre peak to the W of the galaxy, coincident with a knot 
in the optical, is not associated with any noisy bolometers. 

NGC 6120. The submillimetre emission to the W of the 
galaxy, and at the S of the map, is associated with noisy bolometers. 

NGC 6126. The submillimetre source (which we measured as 
a point source) is offset to the S of the optical extent of the galaxy. 
We note that, at minimum contrast, a small satellite/companion ob- 
ject can be seen in this region in the DSS and 2MASS images. 
The 3a peak to the S of the map is not associated with any noisy 
bolometers and is coincident with a small object visible in the op- 
tical. This observation, however, was carried out in poor weather. 

NGC 6131. The submillimetre emission to the very NW of 
the galaxy (beyond the main optical extent) may be associated with 
a noisy bolometer. 

NGC 6137. We estimate that ~ 20% of the 850 ^m flux from 
this elliptical galaxy could be due to synchrotron contamination 
from a radio source associated with the galaxy (Vlahakis et al., in 
prep.). Although only the submillimetre emission to the W of the 
galaxy coincides with a noisy bolometer we note that this observa- 
tion had a tilted sky. 

NGC 6146. We estimate that as much as 80% of the 
850 |im flux from this elliptical may be due to synchrotron contam- 
ination from a radio source associated with the galaxy (Vlahakis et 
al., in prep.). 

NGC 6155 The submillimetre map shows extended emission 
to the S and SE of the galaxy at 850 ]im, coincident with a num- 
ber of small galaxies/condensations in the optical. Marquez et al. 
(1999) find one of the spiral arms in this galaxy is directed towards 
the SE. None of the submillimetre peaks in this map are associated 
with noisy bolometers. 

A large aperture was used to measure all the flux associ- 
ated with this object, and these results are listed in Table 1. How- 
ever at 450 n-m any flux appears confined to the main optical ex- 
tent (though the map at 450 |im is very noisy), and thus the flux 
measurement at 450 iJ.m was made using a smaller aperture. Us- 
ing these values of the 850 iJ-m and 450 [im flux we found a two- 
component SED could not be fitted (Section 3.2); the Saso/Ssso 
ratio is simply too low, most likely because we have measured ex- 
tended emission at 850 nm. Thus we also measured the Ss5o us- 
ing a smaller aperture the same size as used at 450 iJ.m, and find 
Ssso = 0.069 ± 0.013 Jy. For this smaller aperture we find that a 
two-component model can just be fitted to the data, and those pa- 
rameters we list in Table 3. 

NGC 6166 is an elliptical and is located in a very busy field 
- it is the dominant galaxy in the cluster Abell 2199. The presence 
of dust lanes is well documented in the literature. We note that our 
SC ANPI measurements are in good agreement with those of Knapp 
et al. (1989). Using all available radio fluxes from the literature we 
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(a) UGC 148: T=(34,18), n=37 



(b) NGC 99: T=(47,17), n=542 



(c) NGC 803: T=(33,18), n=92 (shown) or 
T=(60,19),n=2597 




(d) NGC 3689: T=(59,23), n=910 



(e) PGC 35952: T=(58,18), n=1859 



(f) NGC 3987: T=(44,22), n=279 






(g) NGC 5953/4; T=(54,21), n=277 



(h) NGC 5980: T=(43,18), n=321 



(i) NGC 6090: T=(55,22), n=122 






(j) NGC 6120: T=(45,24), n=76 



(k) NGC 6155: T=(30,20), n=7 



(1) NGC 6190: T=(56,18), n=2684 



Figure 2. Best-fitting two-component SEDs assuming j3 = 2, fitted to the 60, 100, 450 and 850 (im iiuxes (with the exception of (a), see Section 3.1). Solid 
lines represent the composite two-component SED and dot-dash Unes indicate the warm and cold components. Any additional 170 |jm (ISO) fluxes from the 
Uterature (Stickel et al. 2000, 2004) are also plotted, though (with the exception of (a)) not fitted. Note, captions show the fitted parameters as Usted in Table 3, 
so may be the averaged values (see Section 3.2). 
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(m) IC 5090: T=(52,21), n=346 



(n) IC 1368: T=(55,23), n=110 



(o) NGC 7081: T=(32,20), n=6 






(p) NGC 7442: T=(54,20), n=665 



(q) UGC 12519: TK28,17), n=12 
Figure 2 - continued 



(r) NGC 7722: T=(54,20), n=1224 



estimate that as little as 4% or as much as 100% of the 850 ^.m flux 
from this elliptical may be due to synchrotron contamination from a 
radio source associated with the galaxy (depending whether a spec- 
tral index is assumed constant over the whole galaxy or whether it 
is assumed to have a flatter core) - this is a preliminary analysis 
and discussion of this and the other five ellipticals detected in the 
OS sample will be the subject of a separate paper (Vlahakis et al, 
in prep.)- 

NGC 6173. We measure Sioo=0.20 Jy with SCANPI but the 
detection is unconvincing since the coadds do not agree. Therefore 

we give an upper limit at 100 |j.m in Table 1. 

NGC 6190. Some of the data for this object was very noisy 
and unusable. Consequently the remaining data may not be reliable. 

The submillimetre emission to the W of the galaxy lies in a region 
where there is a noisy bolometer in the 850 iJ.m flux map. Thus the 
apertures used also unavoidably encompass some noisy bolome- 
ters, particularly at 450 |a.m, and results for this object should be 
used with caution. However, the rest of the 850 \im emission in the 
map is not associated with any noisy bolometers, so while the flux 
measurements may be unreliable this does not apply to the emission 
extent, which appears to follow the outer spiral arm structure. 

NGC 7081. The submillimetre emission to the E and W of this 
galaxy are not associated with any noisy bolometers. The emission 
to the N and SE is coincident with regions of the map which are 
only slightly noisy, and since this observation was carried out in 
very good weather it is unlikely that any residual features of this 
noise remain in the S/N map (Section 2.7). Though from the optical 
DSS image only the central region of the galaxy (coincident with 
the main submillimetre peak) is clearly visible, there is evidence 
that this spiral has a very faint extended spiral arm structure. This 



is confirmed by optical images from SuperCOSMOS which clearly 
show very knotty and irregular faint spiral arms coincident with the 
peaks of submillimetre emission to the N and W of the galaxy. 

NGC 7442. The 3a submiUimetre peak to NW of main opti- 
cal extent is coincident with faint optical knots and (unlike the la 
peaks elsewhere in the submillimetre map) is not associated with 
any noisy bolometers. 

NGC 7463. This galaxy is part of a triple system with NGC 
7464 (to the S of NGC 7463) and NGC 7465 (not shown in Fig- 
ure 1). At 850 iim we clearly detect emission from both NGC 7463 
and NGC 7464, which seem to be joined by a bridge of submil- 
limetre emission. The flux listed in Table 1 is for NGC 7463 alone, 
measured in an aperture corresponding to its main optical extent. 
Unfortunately a very noisy bolometer to the SE prevents us measur- 
ing the flux from the eastern half of NGC 7464, but excluding this 
region we measure a flux for the pair of 0.051±0.012 Jy, though 
obviously a lower limit. 

An IRAS source is associated with NGC 7465, which is re- 
solved from the other members of the system at 60 iJ.m (HIRES; 
Aumann, Fowler & Melnyk 1990). Dust properties of this system 
(using SCUBA data observed as part of SLUGS) are studied in de- 
tail by Thomas et al. (2002). 

UGC 12519. The 850 ^.m emission to the NE of this galaxy is 
coincident with a number of small objects seen in the optical and 
is not associated with any noisy bolometers. Although UGC 12519 
is also detected at 450 |i.m the slightly smaller field of view of the 
short array means these NE objects lie just outside the 450 \im map. 

NGC 7722. Along with a very high ratio of cold-to-warm dust 
this object also has a very prominent dust lane, extending over most 
of the NE 'half' of the galaxy. The 850 ia.m emission clearly follows 
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the dust lane evident in the optical. The 2a submillimetre peak to 
the NW of the galaxy is not associated with any noisy bolometers. 



3.2 Spectral fits 

In this section we describe the dust models we fit to the spectral 
energy distributions (SEDs) of the OS sample galaxies and present 
the results of these fits. Comparison of the results of the OS and 
IRS samples is discussed in Section 4. 

DOO found that for the IRS sample the 60 ^.m, 100 ^.m and 
850 |j.m fluxes could be fitted by a single-temperature dust model. 
However, with the addition of the 450 |xm data (DEOl) they found 
that a single dust emission temperature no longer gives an ade- 
quate fit to the data, and that in fact two temperature components 
are needed, in line with the paradigm that there are two main dust 
components in galaxies (Section 1.2). For the OS galaxies we have 
fitted a two-component dust model where there is a 450 |im flux 
available. Since only 19 of the galaxies have 450 iim flux densi- 
ties we have also fitted an isothermal model to the data for all the 
galaxies in the OS sample which were detected at 850 ]xm. 

We fitted two-component dust SEDs to the 60|i.m, 
100 ^m (IRAS), 450 ^m & 850 ^m (SCUBA) fluxes, by minimis- 
ing the sum of the residuals. This two-component model ex- 
presses the emission at a particular frequency as the sum of two 
modified Planck functions ('grey-bodies'), each having a different 
characteristic temperature, such that 

S, = N^x v^B{v, T^) + NcX v^B{v, n) (2) 

for the optically thin regime. Here Ntj, and Nc represent the relative 
masses in the warm and cold components, and Tc the tempera- 
tures, -B(i^, T) the Planck function, and P the dust emissivity index. 
DEOl used the high value for the ratio of S450 /Ssso and the tight 
correlation between Seo /S450 and ^eo / Saso for the IRS galaxies to 
argue that /3 « 2. The OS galaxies follow a similar tight correlation 
(Section 4.1). For the OS sample we find the mean S450 /S85()=8.6 
with (7=3.3, which is slightly higher than found for the IRS sample 
(where Sa5o / Sszo=l .9 with (7=1.6) and with a slightly less tight 
distribution (though still consistent with being produced by the un- 
certainties in the fluxes). Both the OS and IRS values are somewhat 
higher than that found for the Stevens et al. (2005) sample of 14 
local spiral galaxies, where the mean 5450 /5860=5.9 and (7=1.0. 
Since the OS galaxies have a similar high value for this ratio to the 
IRS galaxies we follow DEOl in assuming (3=1. 

We constrained by the IRAS 25 \im flux (the fit was not 
allowed to exceed this value), though we did not actually fit this 
data point, and we allowed To to take any value lower than Tu, . 
This method is the same as that used in DEOl for the IRS sam- 
ple, but while many of the IRS galaxies with 450 |j.m data also had 
fluxes at several other wavelengths in the literature we note that 
for the OS sample galaxies we have only four data points to fit. 
Since this is not enough data points to provide a well-constrained 
fit the values of Xmin may be unrealistically low. In Table 3 we list 
the parameters producing the best fits or, where more than one set 
of parameters produces an acceptable fit, we list an average of all 
those parameters. In practice we find that it is only (and hence 
Nc/Nu,) for which there is sometimes a fairly large range of ac- 
ceptable values, and this is likely due to our not fitting any data 
points below 60 |im. We show all our fitted two-component SEDs 
in Figure 2 (for the best-fitting, not averaged, parameters). Any ad- 
ditional 170 |j.m (ISO) fluxes available from the Uterature (Stickel 
et al. 2000, 2004) are also plotted in Figure 2, though (wifli flie ex- 
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Figure 3. Two representative isothermal SEDs. 



ception of UGC 148 (see Section 3.1)) not fitted; where there are 
several 170 ]im measurements available we plot the mean value. 

We find a mean warm component temperature 
Ttu — 47.4 ± 2.4 K and a mean cold component temperature 
Tc = 20.2 ± 0.5 K. The fitted warm component temperatures are 
in the range 28 < < 59 K and cold component temperatures 
are in the range 17 < Tc < 24 K. Thus, the cold component 
temperature is close to that expected for dust heated by the general 
ISRF (Cox et al. 1986), one of the two components in the current 
paradigm (Section 1.2). 

We find a mean Nc/N^ = 532 ±172 (or higher if we include 
the higher value found for NGC 803 (see notes to Table 3 and Sec- 
tion 3.1)). For the IRS sample DEOl found a large variation in the 
relative contribution of the cold component to the SEDs (described 
by the parameter Nc/Nw); for the OS sample we find an even larger 
variation. Objects NGC 6190 and PGC 35952 in Figure 2, for ex- 
ample, clearly exhibit very 'cold' SEDs with a strikingly prominent 
cold component (with fs 2000 times as much cold dust as warm 
dust). Comparison of the two samples is discussed in detail in Sec- 
tion 4. 

Since 450 |im fluxes are only available for only ^ one third of 
the sample we have in addition fitted single-component SEDs for 
all sources in the OS sample. In these fits we have allowed f3 to vary 
as well as Td, as it is rarely possible to get an acceptable fit with 
8=2. The best fitting Td and 3 are listed in Table 1. We include fit- 
ted parameters only for those objects with detections in all 3 wave- 
bands (60 |im, 100 n-m and 850 |im). The sample mean and error 
in the mean for the best-fitting temperature isT'd = 31.6±0.6K 
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Table 4. Luminosities and masses. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


Name 


log Leo 


log Lsso 


log Lfir 


log Md 


log Mhi 


log Lb 




(WHz~isr~i) 


(WHz~isr~i) 


(Lq) 


(Mq) 


{Mq) 


(Lq) 


UGC 148 


22.80 


21.20 


10.22 


7.05 


9.82 


10.39 


NGC99 


22.57 


21.46 


9.94 


7.17 


10.29 


10.37 


PGC 3563 


22.24 


21.13 


9.76 


6.99 




10.03 


NGC 786 


22.56 


21.34 


9.99 


7.14 




9.94 


NGC 803 


21.69 


20.82 


9.37 


6.76 


9.78 


10.13 


UGC 5129 


21.86 


<20.96 


<9.46 


<7.09 


9.34 


10.01 


NGC 2954 


<21.60 


<20.81 


<9.23 




8.09 


10.25 


UGC 5342 


22.46 


21.03 


9.84 


6.81 




10.40 


PGC 29536 


<22.40 


<21.76 


<9.94 






10.74 


NGC 3209 


<21.99 


<21.13 


<9.66 






10.45 


NGC 3270 


22.58 


21.58 


10.23 


7.53 


10.49 


10.78 


NGC 3323 


22.81 


21.48 


10.23 


7.30 


9.68 


10.12 


NGC 3689 


^22.54 


21.09 


10.09 


7.04 


9.14 


10.29 


UGC 6496 












10.10 


PGC 35952 


22.08 


21.11 


9.61 


6.96 


9.73 


10.07 


NGC 3799/3800^ 


22.93 


21.38 


10.39 


7.27 


9.34 




NGC 3812 


<21.68 


<20.91 


<9.17 






9.95 


NGC 3815 


22.19 


20.96 


9.70 


6.83 


9.62 


10.15 


NGC 3920 


22.21 


20.86 


9.63 


6.68 




9.87 


NGC 3987 


23.20 


21.79 


10.74 


7.72 


9.75 


10.63 


NGC 3997 


22.63 


<20.93 


<9.96 


<7.06 


9.83 


10.30 


NGC 4005 




<20.69 




<6.82 


9.22 


10.35 


NGC 4015 


21.88 


<21.18 


<9.49 


<7.31 






UGC 7115 


<22.17 


21.59 


<9.80 


t7.71 




10.45 


UGC 7157 


<22.15 


<21.27 


<9.65 






10.32 


IC797 


21.72 


20.78 


9.27 


6.63 


8.50 


9.77 


IC800 


21.52 


20.82 


9.07 


6.63 


8.51 


9.67 


NGC 4712 


22.18 


21.50 


9.87 


7.43 


10.18 


10.50 


PGC 47122 


<21.95 


<21.46 


<9.71 


<7.58 




10.32 


MRK 1365 


23.32 


21.20 


10.61 


7.00 


9.23 


10.00 


UGC 8872 


<22.04 


<21.02 


<9.44 






10.29 


UGC 8883 


22.36 


<21.31 


<9.86 


<7.44 




10.00 


UGC 8902 


23.07 


21.81 


10.57 


7.69 




10.80 


IC 979 


**22.27 


21.75 


*9.85 


*7.56 




10.64 


UGC 9110 


U 


<21.21 






9.72 


10.27 


NGC 5522 


22.84 


21.39 


10.22 


7.17 


9.77 


10.51 


NGC 5953/4P 


22.80 


21.22 


10.17 


7.03 


9.32 




NGC 5980 


22.97 


21.84 


10.44 


7.65 




10.53 


IC 1174 


<21.80 


20.95 


<9.19 


t7.08 




10.18 


UGC 10200 


21.95 


<20.10 


<9.21 


<6.22 


9.54 


9.05 


UGC 10205 


22.44 


21.61 


10.10 


7.53 


9.57 


10.55 


NGC 6090 


23.93 


22.06 


11.20 


7.78 


8.82 


10.73 


NGC 6103 


22.97 


21.88 


10.45 


7.71 




10.83 


NGC 6104 


22.77 


<21.59 


<10.35 


<7.71 




10.62 


IC 1211 


<21.78 


21.16 


<9.52 


t7.29 




10.37 


UGC 10325 


22.92 


21.34 


10.35 


7.20 


9.92 




NGC 6127 


<21.57 


21.51 


<9.17 


t7.64 




10.66 


NGC 6120 


23.74 


21.95 


11.11 


7.80 




10.73 


NGC 6126 


<22.37 


21.56 


<9.89 


t7.69 




10.61 


NGC 6131 


22.49 


21.36 


10.07 


7.27 


9.83 


10.37 


NGC 6137 


<22.41 


21.62 


<9.94 


t7.75 




11.04 


NGC 6146 


<22.18 


21.55 


<9.83 


t7.68 




11.01 


NGC 6154 


<21.95 


<21.37 


<9.44 




9.86 


10.38 


NGC 6155 


22.25 


21.04 


9.78 


6.93 


8.95 


9.82 


UGC 10407 


23.28 


21.48 


10.63 


7.32 




10.62 


NGC 6166 


*22.14 


22.00 


10.03 


7.96 




11.30 


NGC 6173 


<22.33 


<21.48 


<9.63 






11.14 


NGC 6189 


22.59 


21.57 


10.19 


7.48 


10.07 


10.68 


NGC 6190 


22.02 


21.24 


9.69 


7.18 


9.48 


9.97 
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Table 4 - continued 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Name 


log Leo 


log Lgso 


log Lfir 


log 


log Mhi 


log Lb 




(WHz~isr~i) 


(WHz~isr~'^) 








(Lq) 


NGC 6185 


22.47 


<21.72 


< 10.06 


<7.85 




10.96 


UGC 10486 


<22.06 


<21.24 


<9.63 






10.31 


NGC 6196 


<22.24 


<21.53 


<9.86 






10.85 


UGC 10500 


**21.85 


<21.09 


<9.56 


<7.22 




10.20 


IC 5090 


23.64 


22.23 


11.09 


8.08 




10.46 


IC 1368 


23.00 


21.07 


10.30 


6.83 




10.14 


NGC 7047 


22.35 


21.45 


9.98 


7.38 


9.08 


10.50 


NGC 7081 


22.49 


20.88 


9.89 


6.72 


9.51 


9.75 


NGC 7280 


<20.80 


<20.34 


<8.51 




8.16 


9.85 


NGC 7442 


22.83 


21.60 


10.32 


7.47 


9.75 


10.48 


NGC 7448 


22.84 


21.19 


10.19 


7.04 


9.75 


10.39 


NGC 7461 


<21.70 


<20.81 


<9.33 






9.90 


NGC 7463 


U 


20.60 




^6.73 


9.33 


10.12 


III ZW 093 


23.26 


<21.99 


< 11.06 


<8.12 


9.95 


11.60 


III ZW 095 


<21.92 


<21.24 


<9.93 






9.90 


UGC 12519 


22.37 


21.36 


9.95 


7.26 


9.53 


10.27 


NGC 7653 


22.59 


21.52 


10.17 


7.43 




10.49 


NGC 7691 


22.15 


<20.82 


<9.68 


<6.95 


9.59 


10.24 


NGC 7711 


<21.60 


<20.86 


<9.20 






10.57 


NGC 7722 


22.31 


21.20 


9.94 


7.15 


9.51 


10.48 



(1) Most commonly used name. 

(2) 60 pra luminosity. 

(3) 850 pm luminosity. 

(4) FIR luminosity, calculated by integrating measured SED from 40-1000 ]Am. 

(5) Dust mass, calculated using a single temperature, T^, as listed in Table 1 (T^ derived from fitted SED 
to the 60, 100 and 850 \mi data points). Upper limits are calculated using T^=2Q K. 

(6) HI mass; refs.: Chamaraux, Balkowski & Fontanelli (1987), Haynes & Giovanelli (1988, 1991), 
Huchtmeier & Richter (1989), Giovanelli & Haynes (1993), Lu et al. (1993), Freudling (1995), DuPrie & 
Schneider (1996), Huchtmeier (1997), Theureau et al. (1998), Haynes et al. (1999). 

(7) Blue luminosity, calculated from corrected blue magnitudes taken from the LEDA database. 

P A close or interacting pair which was resolved by SCUBA; parameters given refer to the combined 

system, as in Table 2. 

* Values should be used with caution (see * notes to Table 1). 

t Object was only detected at 850 yura (and not in either IRAS band), so these dust masses should be used 
with caution; these are all early types; dust masses calculated using Ta=20 K. 
^ Dust mass calculated using T^=2Q K, since no fitted value of T^. 
U Unresolved by IRAS. 

Notes on HI fluxes:- 

NGC 7463 Giovanelli & Haynes (1993) note confused HI profile, many neighbours. 
UGC 12519 HI flux from Giovanelli & Haynes (1993) gives HI mass of log M0=9.65. 
NGC 7691 Haynes et al. (1999) gives log Mq=9.88. 

NGC 5953/4 HI flux from FreudUng (1995); Huchtmeier & Richter (1989) give log M© in the range 

8.82-9.20. 

NGC 3799/3800 flux for NGC 3800 but sources may be confused (Lu et al. 1993). 

NGC 803 Giovanelli & Haynes (1993) give log Mo=9.68, and note optical disk larger than beam. 

NGC 6090 Huchtmeier & Richter (1989) give values up to log Mq=10.24. 

NGC 6131 confused with neighbour. 

NGC 6189 Haynes et al. (1999) gives log M0=1O.26. 

NGC 7081 Huchtmeier & Richter (1989) give values up to log Mq=9.69. 

NGC 4712 Huchtmeier & Richter (1989) give log M0=9.91. 



and for the dust emissivity index ^ = 1.11 ± 0.05. Figure 3 shows 
two representative isothermal SEDs. As an example of the potential 
dangers of fitting single-temperature SEDs, we note one of these 
objects (NGC 99) is also fitted with the two-component model 
(shown in Figure 2). NGC 99 can clearly be well-fitted by both 
an isothermal dust model with very flat p (/3=0.4 in this case) and 



a two-component model with much steeper /3 (/3=2); this is also 
the case for NGC 6190 and PGC 35952 described above. We note 
that the low values of /3 found from the isothermal fits are not the 
true values of f3 but rather is evidence that galaxies across all Hub- 
ble types contain a significant proportion of dust that is colder than 
these fitted temperatures, and it is likely that these objects (as for 
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Figure 4. Colour-colour plot: Seo/Sioo versus Seo/S'gso colours for the optically-selected (this work) and /SA5-selected (DOO) SLUGS (filled and open 
points respectively). 



NGC 99) require a two-component model to adequately describe 
their SED. 



For those galaxies for which we have made two-component 
fits (Table 3) we also calculate the two-component dust mass 
(Md2), using 



3.3 Dust masses 

Dust masses for the OS galaxies are calculated using the measured 
850 |xm fluxes and dust temperatures (Td) from the isothermal fits 
(listed in Table 1) using 



SasoD 



Kd(v)B{v,Td) 



(3) 



where Kd is the dust mass opacity coefficient at 850 i^m, B{v, Td) 
is the Planck function at 850 (im for the temperature Td and D is 
the distance. 

As discussed in DOO we assume a value for Kd(v) of 
0.077m^kg~^, which is consistent with the value derived by James 
et al. (2002) from the global properties of galaxies. Though the true 
value of Kd (i') is uncertain, as long as dust has similar properties 
in all galaxies then our relative dust masses will be correct. The un- 
certainties in the relative dust masses then depend only on errors in 
Ssso and Td. 

Values for dust masses (calculated using Td from our isother- 
mal fits) are given in Table 4. We find a mean dust mass 
Md = (2.34 ± 0.36) X lO'^ M© (where the ± error is the error on 
the mean), which is comparable to that found for the IRS sample 
(DOO). This, together with the fact that for the OS sample we find 
significantly lower values of /3, poses a mmiber of issues. As shown 
by DEOl, if more than one temperature component is present our 
use of a single-temperature model will have given us values of 
P which are lower than actually true, biased our Td estimates to 
higher temperatures, and lead to underestimates of the dust masses. 



Md2 = X 



B{v,Tc) B{v,T^) 



(4) 



where parameters are the same as in Equation 3 and Tc, Tu-, Nc and 
are the fitted two-component parameters as in Equation 2 (and 
listed in Table 3). The mean two-component dust mass is found to 
be Md2 = (4.89 ± 1.20) x lO'^ M©, and the two-component dust 
masses are typically a factor of 2 higher than found from fitting 
single-temperature SEDs, though in some cases (such as NGC 99) 
as much as a factor of 4 higher. 

Given the lack of CO measurements for the OS sample galax- 
ies, one potential problem with the above estimates of dust mass 
is any contribution to the SCUBA 850 \im measurements by CO(3- 
2) line emission. Seaquist et al. (2004) find, for a representative 
subsample of the IRS SLUGS galaxies from DOO, that contami- 
nation of 850 n-m SCUBA fluxes by CO(3-2) reduces the average 
dust mass by 25-38%, though this does not affect the shape of the 
dust mass function derived using the IRS SLUGS sample in DOO. 
However, the OS galaxies are relatively faint submillimetre sources 
compared with the IRS sample. From the fractional contribution 
of CO(3-2) line emission derived by Seaquist et al. (a linear flt to 
the plot of SCUBA-equivalent flux produced by the CO line ver- 
sus SCUBA flux) we estimate that for the OS sample the CO line 
contribution to the 850 \im flux is small and is well within the un- 
certainties on the 850 \im fluxes we give in Table 1. 
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Figure 5. Colour-colour plot: Seo / S450 versus Sao / S'gso colours for the optically-selected (this work) and /R45-selected (DOO) SLUGS (filled and open 
points respectively). 



3.4 Gas masses 

The neutral hydrogen masses listed in Table 4 were calculated from 
HI fluxes taken from the literature^ using 

Mhi = 2.356 x W^D^Shi (5) 

where D is in Mpc and Shi is in Jy km s~^. 

Only a small handful of objects in the OS sample had CO 
fluxes in the literature, and so in this work we will not present any 
molecular gas masses. 

3.5 Far-infrared luminosities 

The FIR luminosity (Lfir) is usually calculated using 

FIR = 1.26 X 10""(2.58S60 + Swo) 

and 

L/„. = AttD'-^ X FIB, X C (6) 

as described in the Appendix of Catalogued Galaxies and Quasars 
Observed in the IRAS Survey (Version 2, 1989), where Seo and 
S'loo are the 60 |im and 100 ^.m IRAS fluxes, D is the distance, and 
C is a colour-correction factor dependant on the ratio 5*60 / Sioo and 
the assumed emissivity index. The purpose of this correction factor 
is to account for emission outside the IRAS bands, and is explained 
by Helou et al. (1988). 

However, since we have submillimetre fluxes we can use our 
derived Td and /3 to integrate the total flux under the SED out to 
1000 |im. This method gives more accurate values of Lfir since it 

^ See notes to Table 4. 



makes no general assumptions. We list in Table 4 Lfir calculated 
using this method and our fitted isothermal SEDs; Lfir values cal- 
culated using our two-component SEDs are Usted in Table 3. 

3.6 Optical luminosities 

The blue luminosities given in Table 4 are converted (using 
Mb0=5.48) from blue apparent magnitudes taken from the Lyon- 
Meudon Extragalactic Database (LEDA; Paturel et al. 1989, 2003) 
which have already been corrected for galactic extinction, internal 
extinction and k-correction. 

4 THE SUBMILLIMETRE PROPERTIES OF GALAXIES 

4.1 Optical selection versus IR selection 

Figures 4 and 5 show the OS and IRS galaxies plotted on two- 
colour diagrams (filled and open symbols respectively). The IRS 
and OS galaxies clearly have different distributions, and in particu- 
lar there are OS galaxies in parts of the diagram where there are no 
IRS galaxies. In Figure 4 ~ 50% of the OS galaxies are in a region 
of the colour-colour diagram completely unoccupied by IRS galax- 
ies. This shows there are galaxies 'missing' from IR samples, with 
important implications for the submillimetre LF (Section 5). 

Figure 5 shows the S6o/>S'450 versus Sw/ SsrA) colour- 
colour plot for the OS sample objects and IRS sample objects 
which have 450 ^m fluxes. We confirm the very tight correla- 
tion found by DEOl (here the correlation coefficient rs = 0.96, 
significance = 9.20e-21)), and the scatter for the OS sample may 
be completely explained by the uncertainties on the fluxes. Impor- 
tantly, this relationship holds for all the objects in the OS sample 
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(a) (b) 
Figure 6. Distributions of (a) /? values and (b) values for the optically- and fflA5-selected SLUGS (line-fiUed and shaded histograms respectively). 



for which we have 450 yita fluxes, which include a wide range of 
galaxy types (t-type=0 to 10) and with Lfir ranging over 2 orders 
of magnitude. The (least-squares) best-fitting line to the combined 
OS -I- IRS samples shown in Figure 5 is given by 

log(S6o/S46o) = (1-03 ± 0.05) log(S6o/S85o) - (0.955 ± 0.070) 

(or re-written 5'6o/5'450 = 0.119(5'6o/S'85o)^ °^) and is very 
similar to that found by DEOl, confirming the finding for the 
IRS sample that, within the uncertainties, the ratio 5450/5*850 is 
constant. DEOl conclude, from the results of simulations of the 
450/850 i^m flux ratio and from the fitted (3 values for those galax- 
ies whose SEDs require a cold component, that 13 ^ 2 for all galax- 
ies, and that therefore the cold dust component in all galaxies has 
a similar temperature (Tc ~ 20-21 K). The fact that we also find 
the 5460 /5860 ratio constant for the OS sample suggests that these 
conclusions are true for all Hubble types (only t-types<0 are un- 
represented in the OS sub-sample with 450 |j.m data). 

The positions of the OS galaxies in the colour diagrams sug- 
gest there is more cold dust in the OS galaxies than in the IRS 
galaxies. We can investigate this further with the results of our spec- 
tral fits. Figure 6(a) shows the comparison between the distribution 
of 13 values (found from the isothermal fits) for the OS and IRS 
samples. We find OS sample galaxies with (3 values lower than any 
found in the IRS sample. The two-sided Kolmogorov-Smimov (K- 
S) test shows that the distributions of the two samples are signif- 
icantly different (the probability that the two samples come from 
the same distribution function is only 1.8e-5). Though this clearly 
demonstrates that the properties of the dust in the OS and IRS sam- 
ples are different, rather than interpreting this as a physical differ- 
ence in the emissivity behaviour of the grains (B) we believe that it 
is a difference in the two samples' ratios of cold/warm dust. 

Figure 6(b) shows the comparison between the distribution of 
dust temperatures (from isothermal fits) for the OS and IRS sam- 
ples. We note that the OS sample has consistently colder Td com- 
pared to the IRS sample. Once again using a K-S test we find that 
the OS and IRS sample dust temperatures do not have the same dis- 
tribution, with the probability of the two samples coming from the 
same distribution function being only 1.41e-4. 



For those objects in the OS and IRS samples for which two- 
component fits were possible, the distributions of the warm and 
cold component temperatures (T^, and Tc) for the two samples 
are shown in Figures 7(a) and 7(b) respectively. The distribu- 
tions of Tc for the OS and IRS samples are statistically indistin- 
guishable, while conversely the distributions of T-u, are not simi- 
lar (probability of same distribution is 0.03). While the mean cold 
component temperature for the OS sample (Tc = 20.2 ± 0.5 K) 
is very similar to the value found for the IRS sample (mean 
Tc = 20.1 ± 0.4 K), the mean warm component temperature is 
rather higher (TU = 47.4 ± 2.4 K for the OS sample as opposed 
to = 39.3 ± 1.4 K for the IRS sample). 

Figure 8 shows the distribution of Nc/N^ for the OS and 
IRS samples. The OS and IRS samples clearly have different dis- 
tributions - the (K-S test) probability of the two samples hav- 
ing the same distribution is 8.4e-4. For the OS sample the mean 
Nc/N.^ = 532 ± 172 (or higher, see Section 3.2), for the IRS sam- 
ple the mean Nc/N^ = 38 ± 11. For the OS sample there is a 
much larger range of Nc/Nw than for the IRS sample. Interest- 
ingly, few of the OS objects even have a NJN^n low enough to 
fall within the range found for the IRS sample, strongly suggesting 
a prevalence of cold dust in the OS sample compared to the IRS 
sample. 

The large difference between the distributions of Nc/N^ for 
the OS and IRS galaxies implies that most OS sample galaxies con- 
tain much larger proportions of cold dust relative to warm dust than 
found for the IRS galaxies, additional evidence that IRAS missed a 
population of cold-dust-dominated objects. The similarity of the 
temperature of the cold component for the OS and IRS sample 
and the difference in the distribution of N^/N-u) supports the cur- 
rent paradigm for dust in galaxies. An alternative model for dust in 
galaxies would be one in which IRAS galaxies are ones in which the 
general ISRF is more intense, and therefore the majority of dust is 
hotter. The similarity of Tc for the different samples argues against 
this and suggests that most dust in all galaxies is relatively cold 
and has a similar temperature. The temperature differences between 
galaxies arise from a second dust component, presumably the dust 
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Figure 8. Distribution of log{Nc/Niu) for tlie OS and IRS samples (line- 
filled and shaded histograms respectively). 
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Figure 9. Nc/Nw versus 60 ^mi luminosity for the OS and IRS samples 
(filled and open points respectively). 



Figure 7. Distributions of warm component (a) and cold component (b) 
temperatures for the OS and IRS samples (line-filled and shaded histograms 
respectively). 



in regions of intense star formation. Our results for the OS sample 
indicate that the ratio of the mass of dust in this second component 
to the mass of dust in the first component can vary by roughly a 
factor of 1000. There are two other pieces of evidence in favour of 
the two-component model. First, the ISO 170 iJ-m fiux densities that 
exist for 3 of our two-component-fitted galaxies (Stickel et al. 2004; 
Section 3.2) agree very well with our model SEDs (we did not use 
these data in making our fits, with one exception; see Section 3.2). 
Second, the ratio of the mass of cold dust to the mass of warm dust 
correlates inversely with 60 |j.m luminosity (Figure 9; Ts = —0.41, 
significance = 1.24e-2); in the two-component model one might ex- 
pect the most luminous IRAS sources to be dominated by the warm 
component. 

The difference in the distributions of Tw does not, however, fit 
in with this general picture. In the two-component model one would 
expect Tyj and Tc to be constants, with the only thing changing be- 
tween galaxies being the proportion of cold and warm dust. The 



difference in the distributions of T^, may indicate that this model 
is too simplistic. Two things may be relevant here. First, as can be 
seen in Figure 2 it is those OS galaxies with very prominent cold 
components which typically account for the highest warm compo- 
nent temperatures (for example PGC 35952 or NGC 6090). Second, 
the model SEDs with high values of also generally provide a 
good fit to the 25 ^m flux density, whereas the model values with 
low values of tend to underestimate the 25 |i.m flux density. This 
last point suggests that to fully understand dust in galaxies one can- 
not ignore the measurements at wavelengths < 60 [im; however, if 
we did include these measurements we would then definitely need 
more than two dust components. This is clearly demonstrated by 
Sievers et al. (1994) who, for NGC 3627, fit a three-component 
model. A two-component model is nonetheless adequate for our 
purposes, since we are interested in the cold component rather than 
a third hot component. 



4.2 Submillimetre properties along the Hubble sequence 

In this section we investigate the submillimetre properties of galax- 
ies as a function of Hubble type (f). We first compare the distri- 
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Figure 10. Distribution of Hubble types for the OS sample 
850 (jm detections (upper panel) and non-detections (lower panel). 



butions of Hubble type for the detections (D) and non-detections 
(ND) in our OS sample (Figure 10). We use the K-S test to find 
that the probability of their having the same distribution is ~ 2%. 
Thus Figure 10 suggests that early-type galaxies are less likely to 
be submillimetre sources than later types. 

To investigate this apparent morphological difference further, 
we estimated the submillimetre luminosity distributions of early- 
and late-type galaxies. A major complication is the large number 
of upper limits. We used the Kaplan-Meier estimator (Wall & Jenk- 
ins 2003) to incorporate information from both the upper limits and 
the measurements. We defined early-type galaxies as all those with 
t ^ 1 and late-type galaxies as those with t > 1. We used this 
division, because the greatest difference between the cumulative 
distributions of Hubble type for detected and non-detected galax- 
ies (Figure 10) was found at t=l. Figure 11 shows the cumulative 
luminosity distributions estimated in this way for the early-type 
and late-type galaxies. There appears to be a tendency for the late- 
type galaxies to be more luminous submillimetre sources. However, 
the tendency is not very strong. We also used the ASURV statisti- 
cal package for censored data (Feigelson & Nelson 1985) to com- 
pare the results for the two samples, using the Gehan test and the 
log-rank test (see Wall & Jenkins 2003). We found a marginally 
significant (10%) difference using the log-rank test but no signif- 
icant difference using the Gehan test. Figure 12 shows a plot of 
850 nm luminosity versus optical luminosity. For clarity we sim- 
ply divide our sample into 4 broad groups based on the galaxies' 
t-type parameter given in LEDA (which uses the standard numer- 
ical codes for the de Vaucoulers morphological type, as defined 
in RC2): E-SO (t=-5 to 0), Early-type spirals (t=l to 4), S? (t=5) 
and Late-type spirals (t=6 to 10). The different Hubble types show 
similar relationships. On further inspection of the data, the more 



Figure 11. Cumulative luminosity distributions for early-type galaxies 
(solid line) and late-type galaxies (dot-dashed). The maximum values for 
both samples are less than one because of the upper limi ts that fall below 
the lowest actual measurement. 




Lb / Lo 

Figure 12. 850 |am luminosity versus optical luminosity Lg for the OS 
sample, with different Hubble types indicated by different symbols: E- 
S0(t=-5 to 0), Early-type spirals(t=l to 4), S?(t=5), Late-type spirals(t=6 
to 10): circles, triangles, stars, and squares respectively. The 6 detected el- 
lipticals are highlighted as open circles. 
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Figure 13. Distribution of (a) /3 values and (b) values for the OS SLUGS, with different Hubble types (as given in LEDA according to the RC2 code, and 
listed in Table 1) indicated by different shaded regions: E-SO (t=-5 to 0), Early-type spirals (t=l to 4), S? (t=5), Late-type spirals (t=6 to 10). ((b): Inserted 
panel: mean for the OS SLUGS for different Hubble types, with error bars of error on the mean (bins with only 1 source are not plotted)). 



marked dependence on Hubble type visible in Figvire 10 appears to 
be at least partly caused by the early-type galaxies being observed 
in worse conditions. In summary, there appears to be some differ- 
ence in submillimetre properties as one moves along the Hubble 
sequence, but it is not very strong. 

We can also use the results of our spectral fits to investigate 
whether there are any trends with Hubble type. As above, we sim- 
ply divide our sample into 4 broad groups based on the galaxies' 
t-type. Figures 13(a) and 13(b) show the distributions of /3 and Td 
(derived from our single-component fits) for the OS sample. We 
note that the objects of each type appear fairly evenly distributed 
across the bins from ,5 = to 2 (Figure 13(a)), and in order to test 
this statistically we divide the sample into two broad groups: early 
types (—5 < t-type ^ 4) and late types (5 ^ t-type ^ 10), and 
perform a K-S test on the two groups. We find that the distributions 
of the early and late type groups are not significantly different. The 
distribution of isothermal dust temperatures appears similar for all 
Hubble types (Figure 13(b)); we find no significant differences be- 
tween the early and late tj^jes. We also investigated the distribu- 
tions of the warm and cold component temperatures found from 
our two-component fits to look for any differences between early 
and late types; for example Popescu et al. (2002) find a tendency 
for the temperatures of the cold dust component to become colder 
for later t3T5es. We divided our 18 two-component fitted temper- 
atures into Hubble types as in Popescu et al. (2002) and, due to 
our smaller number of sources, also into two broad groups of early 
(0 ^ t-type ^ 4) and late (6 ^ t-type ^ 10) types, and compared 
the overall distributions and the median Tc for each t3fpe grouping. 
We found no differences between either the overall distributions or 
the median values of Tc or T„ for the early and late types, though 
we note the limitations of such a small sample. 

4.3 Ellipticals 

It was once thought that ellipticals were entirely devoid of dust and 
gas, but optical absorption studies now show that dust is usually 
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Figure 14. Dust mass versus HI mass for the OS and IRS samples (filled 
and open circles respectively). 



present (Goudfrooij et al. 1994; van Dokkum & Franx 1995). Fur- 
thermore, dust masses for the ~ 15% of ellipticals detected by IRAS 
(Bregman et al. 1998) have been found to be as much as a factor of 
10-100 higher when estimated from their FIR emission compared 
to estimates from optical absorption (Goudfrooij & de Jong 1995). 

At 850 H-m we detect 6 ellipticals, from a total of 1 1 ellipti- 
cals in the OS sample, and find them to have dust masses in ex- 
cess of 10^ Mq. However, a literature search revealed that for 4 
of the 6 detections there are radio sources. We have used the ra- 
dio data to estimate the contribution of synchrotron emission at 
850 \im. These estimates are often very uncertain because of the 
limited number of flux measurements available between 1.4GHz 
and 850 ^.m (353GHz). However, in some cases (Section 3.1) it is 
clear that some or aU of the 850 ^.m emission may be synchrotron 
radiation. We are currently investigating ellipticals further with 
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SCUBA observations of a larger sample. This will be the subject 
of a separate paper (Vlahakis et al., in prep.). 

4.4 The relationship between gas and dust 

In DOO we found that both the mass of atomic gas and the mass of 
molecular gas are correlated with dust mass, but the correlation is 
tighter for the molecular gas. There are virtually no CO measure- 
ments for the OS sample, so here we have only estimated the mass 
of atomic gas. We compared the dust mass (Md, calculated using 
dust temperatures from the isothermal fits) to the HI mass for the 
OS sample (Figure 14) and find a very weak correlation. Though 
the correlation for the OS sample alone is very weak it is nonethe- 
less consistent with the correlation found by DOO for the IRS sam- 
ple; most of the OS points lie within the region covered by the IRS 
points, but though they cover the same range of HI masses we note 
that we do not have any HI masses for our OS sample objects with 
the higher dust masses. The weakness of the correlation for the OS 
sample is therefore likely due simply to the small number of OS 
sample 850 i^m detections for which we have HI data (28 objects). 

The mean neutral gas-to-dust ratio for the OS sample is 
Mhi /Md=395±71, where the error given is the error on the mean. 
The (neutral + molecular) gas-to-dust ratios for the IRS SLUGS 
sample and the Devereux & Young (1990; herein DY90) sam- 
ple of spiral galaxies are respectively Mh2+hi /Md=5Sl±43 and 
MH2+Jfi/A#d=1080±70, but since for the OS sample we have 
no CO measurements and therefore no measure of the mass of 
molecular hydrogen we can at this stage only compare the neu- 
tral gas-to-dust ratio for the OS sample. We therefore compare 
our OS value to mean neutral gas-to-dust ratios which we calcu- 
late, for the IRS sample and the DY90 sample respectively, to be 
Mh I /Ma=305±24 and MH//Md=2089±341. There is a large 
difference between the values for both SLUGS samples and the 
value determined by DY90. This is almost certainly due to the fact 
that the DY90 dust masses were estimated from IRAS fluxes and 
therefore, for the reasons described in Section 1, will have 'missed' 
the cold dust. 

There is also a difference between the SLUGS values and the 

Galactic value of 160 for the (neutral + molecular) gas-to-dust ra- 
tio (the value derived from Sodroski et al. (2004) by DOO). The 
neutral gas-to-dust ratios for both the SLUGS samples are at least 
a factor of 2 larger than this Galactic value, and as shown by DOO 
when the molecular gas is included the value of the gas-to-dust ra- 
tio for the IRS sample is more than 3 times larger than the Galactic 
value. DOO attribute this discrepancy to a missed 'cold dust' compo- 
nent (Td ^ 20 K) in the IRS sample. We have already noted in this 
paper that the single-temperature fits lead to dust masses approxi- 
mately a factor of 2 lower than the more realistic two-component 
fits (Section 3.3). Using the dust masses calculated using our two- 
component fits (Md2', Table 3), for the 13 galaxies for which there 
are HI masses we find the mean neutral gas-to-dust ratio for the 
OS sample is then Mhi /Md2=192±44. This is in good agreement 
with the Galactic value, although if there is a significant amount of 
molecular gas this value would obviously be higher. 



5 LUMINOSITY AND DUST MASS FUNCTIONS 

The 'accessible volume' method (Avni & Bahcall 1980) will, in 
principle, produce unbiased estimates of the submillimetre lumi- 
nosity function (LF) and dust mass function (DMF) provided that 
no population of galaxies is unrepresented by the sample used to 



derive the LF and DMF. In Paper I (DOO) we produced a first es- 
timate of the LF and DMF from the IRS sample. However, since 
our new observations of the OS sample have shown the existence 
of a population of galaxies with low values of the Seo /Sioo and 
Sgo/Ssso flux ratios (Figure 4 and discussion in Section 3.2) of 
which there is not a single representative in the IRS sample, our 
earlier estimates of the LF and DMF are likely to be biased. 

In this section we use our new (OS sample) results to produce 
new estimates of the submillimetre LF and DMF. 



5.1 Method 

We derive the local submillimetre LF and DMF by two different 
methods: 1) directly from the OS SLUGS sample, and 2) by ex- 
trapolating the spectral energy distributions of the galaxies in the 
IRAS PSCz catalogue out to 850 li-m. The PSCz catalogue (Saun- 
ders et al. 2000) is a complete redshift survey of ~ 15000 IRAS 
galaxies in the IRAS Point Source Catalogue. Serjeant & Harrison 
(2005; herein SH05) used the PSCz galaxies and the IRS SLUGS 
submm:far-IR two-colour relation to extrapolate the SEDs of the 
PSCz galaxies out to 850 ia.m and produce an 850 [xm LF. Impor- 
tantly, this method allows us to probe a wider range of luminosities 
than probed directly by the SLUGS samples. 
We estimate the LF for both methods using 

$(L)AL = J2y^ (7) 

i 

(Avni & Bahcall 1980). Here $(L)AL is the number density of 
objects (Mpc~^) in the luminosity range L to L + AL, the sum- 
mation is over all the objects in the sample lying within this lu- 
minosity range, and Vf is the accessible volume of the ith object 
in the sample. Throughout we use an Ho of 75 km s~^Mpc~^ and 
a 'concordance' universe with Om=0.3 and f2A=0.7. We estimate 
the dust mass function (the space density of galaxies as a function 
of dust mass) in the same way as the LF, substituting dust mass for 
luminosity in Equation 7. The details of these two methods, here- 
after referred to as 'directly measured' and 'PSCz-extrapolated', 
are discussed in Sections 5.1.1 and 5.1.2 respectively. 



5. 1. 1 Directly measured 850 fim luminosity function and dust 
mass function 

We calculated the directly measured LF and DMF from the 52 ob- 
jects in the OS sample which were detected at 850 |im. For the 
DMF we use the dust masses listed in Table 4, which were cal- 
culated using the isothermal SED-fitted temperatures or, where no 
fit was made (11 objects), using a dust temperature of 20K. For 
the OS sample the accessible volume is the maximum volume in 
which the object would still be detected at 850 ^im and still be in- 
cluded in the CfA sample. Since objects with cz < 1900 km s~^ 
were excluded from our sample this volume is not included in our 
calculation of Vi. When calculating the maximum redshift at which 
an object would still be detected at 850 |xm we used the noise ap- 
propriate for the observation of that object. We corrected the LF by 
the factor 97/81 to account for the CfA galaxies we did not observe 
at all at 850 jim (Section 2.1). 

The corrected directly measured 850 ^.m LF and DMF are 
shown as star symbols in Figures 15 and 16 respectively, and are 
given in tabular form in Table 5. The errors on the directly mea- 
smed LF and DMF are standard Poisson errors. One effect that 
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Table 5. Directly measured OS SLUGS luminosity and dust mass functions 



Table 6. PSCz-extrapolated luminosity function 
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850 |am dust mass function 
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may lead to our estimates of the LF and DMF being slight under- 
estimates is that we noticed that the OS galaxies not detected at 
850 ^tm were generally observed under worse weather conditions 
than the sources that were detected. 



J. 1.2 IRAS PSCz-extrapolated 850 nm luminosity function and 

dust mass function 

In order to better constrain the LF at the lower luminosity end more 
data points are needed, probing a wider range of luminosities than 
probed directly by the SLUGS samples. We acliieve this using a 
method described by SH05, whereby the 850 H-m LF is determined 
by extrapolating the spectral energy distributions of the ~ 15000 
IRAS PSCz survey galaxies (Saunders et al. 2000) out to 850 \im. 
Since for the two SLUGS samples we find a strong correlation be- 
tween the Seo /'S'loo and Seo /S'sso colours (Figure 4) we can use 
a linear fit to this colour-colour relation to make the extrapolation 
from 60 nm to 850 \im fiux density. SH05 derived the submm:far- 
IR two-colour relationship from the IRS SLUGS sample. However, 
we have shown in this paper that the OS and IRS samples have 
quite different properties. In order to determine the sensitivity of 





<P(L) 


—down 

% 


up 


FT 1 1 \ 

Hz^sr^) 


(Mpc-3dex-'^) 


(Mpc 3(Jex 


^) 


18.52 


3.42e-02 


2.42e-02 


2.74e-02 


18.75 


6.30e-02 


2.38e-02 


2.83e-02 


18.99 


3.90e-02 


1.62e-02 


9.45e-03 


19.23 


3.20e-02 


1.06e-02 


6.17e-03 


19.47 


2.35e-02 


3.50e-03 


7.86e-03 


19.70 


3.08e-02 


8.14e-03 


3.42e-03 


19.94 


1.85e-02 


2.81e-03 


5.72e-03 


20.18 


1.26e-02 


1.98e-03 


1.34e-03 


20.42 


1.16e-02 


1.14e-03 


6.74e-04 


20.65 


1.02e-02 


1.70e-03 


4.41e-04 


20.89 


6.67e-03 


1.12e-03 


8.25e-04 


21.13 


4.30e-03 


6.77e-04 


3.01e-04 


21.36 


2.73e-03 


7.59e-04 


1.63e -04 


21.60 


1.34e-03 


4.61e-04 


l.OOe-04 


21.84 


4.43e-04 


1.75e-04 


1.36e-04 


22.08 


1.17e-04 


6.67e-05 


2.36e-05 


22.31 


1.85e-05 


7.86e-06 


1.45e-05 


22.55 


4.27e-06 


2.64e-06 


3.81e-07 


22.79 


2.91e-07 


1.28e-07 


9.39e-07 


23.03 


9.86e-08 


5.62e-08 


3.49e-08 



(WHz-isr-i) (Mpc-3dex-i) 



-1.38 



-t-0.02 
0,03 



3.73 



+0.29 
0,32 



X 10*= 



7+0.41 



xl 



1.0 



Table 7. Linear fit parameters for the SLUGS colom-colour 

(log(S'6o/5'ioo) vs log(56o/5g5o)) shown in Figure 4. 



plot 



SLUGS data fitted 


linear fit (y=mx+c) 




m 


c 


OPT+/R45 


0.365 ± 0.014 


-0.881 ±0.024 


OPT 


0.296 ± 0.031 


-0.797 ± 0.039 


IRAS 


0.421 ± 0.023 


-0.981 ± 0.042 



the LF/DMF to the colour relationship we have derived colour rela- 
tionships for the combined OS + IRS sample, the OS sample alone, 
and the IRS sample alone (Table 7). 

In order to produce unbiased estimates of the LF and DMF 
we have excluded some PSCz galaxies. Firstly we exclude all 
those objects that do not have redshifts, those that have veloci- 
ties < 300 km s"'^ (to ensure that peculiar velocities are imimpor- 
tant), and those that have redshifts > 0.2 (thus excluding any ultra- 
luminous, high-redshift objects). We then exclude all objects with 
upper limits at 100 |j.m, since for these objects we cannot apply the 
SH05 method. Finally, we use the IRAS Point Source Catalogue 
flags, as listed in the PSCz catalogue, to exclude sources which 
are likely to be either solely or strongly contaminated by Galac- 
tic cirrus. It is important to exclude these sources because they are 
very cold sources and so potentially can have a large effect on the 
850 ixm LF. If two or more of the flags indicate Galactic cirrus (us- 
ing flag value Umits indicated in the IRAS Explanatory Supplement) 



32 Catherine Vlahakis, Loretta Dunne, Stephen Bales 



T 




L850 (W Hz-' sr-l) 

Figure 15. PSCz-extrapolated 850 |jm luminosity function (filled circles) with best-fitting Schechter function (soUd Une). The parameters for the Schechter 
function are a = —1.38, L* = 3.7 x 10^^ WHz~^sr~^. Also shown are the directiy measured 850 (Jm luminosity function for the OS SLUGS sample 
(filled stars) with best-fitting Schechter function (dashed line) and the results for the IRS SLUGS sample from Dunne et al. (2000) (open triangles and dotted 
Une). 



we exclude that object. As a check on the validity of this method 
we inspected by eye (using the IRSA ISSA Image Server) a sample 
of ~40 objects randomly chosen from those excluded as Galactic 
cirrus, and a further sample of ~40 objects randomly chosen from 
those that made it into our final sample. We found that 98% of the 
sources with cirrus flags and 7% of the sources without cirrus flags 
showed signs of significant cirrus, although for two thirds of the 
sources with cirrus flags there still appeared to be a genuine source 
present. In total, from the ~ 14500 galaxies with redshifts in the 
IRAS PSCz catalogue we exclude ~4300 objects because of either 
100 [im upper limits or Galactic cirrus. This leaves 10252 galaxies 
in our PSCz-selected sample. 

For the PSCz-extrapolated sample the accessible volume is 
the maximimi volume in which the object could still be seen and 
still be included in the IRAS PSCz catalogue. Since objects with 
cz < 300kms~^ were excluded from our sample this volume is 
not included in our calculation of Vi. For the PSCz-extrapolated 
DMF the dust masses were calculated using Ta derived from the 
IRAS 100/60 colour and 0=2. 

For completeness, the effect of excluding real 60 ^.m sources 
must be taken into account by applying a correction factor to the 
LF and DMF. This correction factor will be uncertain, since some 
excluded sources will be real and some not. Therefore we correct 
using our best estimate of real sources as follows. We corrected for 
two thirds of the sources we excluded as being contaminated by 
cirrus. The correct correction factor for the sources that were ex- 
cluded because they have 100 p.m upper limits is even more uncer- 
tain. These are probably all genuine sources, but they will generally 
have warmer colours than the sources that were not excluded. We 
arbitrarily corrected for 50% of these. Including the correction for 
~100 sources without redshifts, the final correction factor for ex- 
cluded sources is 1.27. This is obviously very tmcertain, however at 



the most it could be 1.43 and at the least it could be 1.00. This pro- 
duces maximum errors of -|-13% and —21% on the LF and DMF 

in addition to the errors described below. We made a correction for 
evolution out to z=0.2 using a density evolution oc (1 -|- zY (Saun- 
ders et al. 1990). We confirmed that the strength assumed for the 
evolution made virtually no difference to our results. 

The PSCz-extrapolated 850 \im LF and DMF are shown as 
filled circles in Figures 15 and 16 respectively, and are given in 
tabular form in Tables 6 and 8. For comparison we also produce 
a 'cold' PSCz-extrapolated DMF, produced as above but with dust 
masses calculated using rd=20K and /3=2; this is shown as open 
circles in Figure 16(a) and Usted in Table 8. 

While the errors on the directly measured LF and DMF are 
standard Poisson errors, the errors on the PSCz-extrapolated LF 
and DMF are derived from a combination of Poisson errors and 
the errors resulting from the fact that the 850 ^.m luminosities have 
been derived using the best-fitting linear relation to our SLUGS 
colour-colour plot (Figure 4). In order to take into account how our 
'choice' of linear fit affects the LF we produce, we additionally 
generate two 'extremes' of the PSCz-extrapolated LF and DMF us- 
ing two alternative fits to the SLUGS colour-colour plot: 1) a fit to 
the OS data only, and 2) a fit to the IRS data only (linear fit param- 
eters listed in Table 7). We then use the maximum difference be- 
tween these 'extreme' LF values and our actual PSCz-extrapolated 
LF data points as the errors on our LF due to our 'choice' of colour- 
colour linear relation. We then also take into account the number 
statistics, and thus add in quadrature the standard Poisson errors 
and the 'choice of colour-colour fit' errors to obtain our total er- 
rors listed in Tables 6 and 8. In addition to these errors there are, at 
most, upper and lower errors of -i-13% and —21% from our choice 
of correction factors. 
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Figure 16. (a) PSCz-extrapolated dustmass function (filled circles) with best-fitting Schechter function (solid line). The dust masses were calculated using 
Trf derived from the IRAS 100/60 colour and 13=2. The parameters for the Schechter function are a = -1.34, M, = 2.7 x lO'^M©. The dashed hne 
and open circles are for a 'cold dustmass function' in which dust masses are calculated using T^=2QK and /3=2. The best-fitting Schechter parameters are 
a = —1.39, M* = 5.3 x IO^Mq. (b) Directly measured dustmass function for the OS SLUGS sample (filled stars) with best-fitting Schechter fimction 
(dashed line). The dust masses were calculated using from isothermal SED fitting. Also shown are the results for the IRS SLUGS sample from Duime et 
al. (2000) (open triangles and dotted Une). The filled circles and sohd Une show the PSCz-extrapolated dustmass function as in (a). 
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Table 8. PSCz-extrapolated dustmass functions 



PSCz-extrapolated single temperature dustmass function 
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5.2 Results and discussion 

The directly measured OS LF and PSCz-extrapolated LF agree 
remarkably well over the range of luminosities covered by the 
SLUGS samples, yet we find that in comparison the IRS sample 
of DOO (plotted as triangles in Figures 15 and 16) consistently un- 
derestimates the submillimetre LF by a factor of 2 and the DMF 
by a factor of 4. The fact that we see this underestimate compared 
to our OS sample, which by definition should be free from any 
dust temperature selection effects, is strong evidence that a popu- 
lation of 'cold' dusty galaxies was indeed 'missed' by IRAS and 
that therefore the IRS sample was missing '^half the galaxies. The 
bigger difference between the DMFs is probably due to the fact 
that, unlike the IRS sample, for the OS sample we do not have fit- 
ted isothermal SEDs for all galaxies and therefore have calculated 
dust masses using an assumed Td=20K for ~20% of the sample 
(Section 5.1.1). 

We fit both the directly measured and PSCz-extrapolated 
850 |im LFs and DMFs with Schechter functions of the form 

L 



(Press & Schechter 1974; Schechter 1975). The best-fitting param- 
eters for the PSCz-extrapolated 850 ^.m LF and DMF are listed in 
Tables 6 and 8 respectively, along with the reduced chi-squared val- 
ues (xi) for the fits; likewise best-fitting parameters for the directly 
measured 850 jim LF and DMF are shown in Table 5. 

We find that both the directly measured and PSCz- 
extrapolated LFs and DMFs are well-fitted by Schechter functions. 
For the PSCz-extrapolated LF and DMF the best-fitting Schechter 
function (a = — 1.38) fits the data points extremely well across 
most of the luminosity range - however, we note that the PSCz- 
extrapolated functions are much less well fitted at the high lumi- 
nosity end. Investigation of the 3 or 4 high end luminosity bins 
has found several anomalies for the objects in these bins, the most 
striking of which is the fact that in each bin there is typically a small 
number of objects with accessible volumes 2 or 3 orders of magni- 
tude lower than the rest of the objects in that bin, and thus it is these 
few objects in each of these bins which are the main contributors to 
the high space density. 

There are many possible explanations for the excess at the 
high luminosity end. One possible explanation could be that the 
objects in these bins are multiple systems. At larger distances IRAS 
galaxies are mostly very luminous starbursts and are frequently 
in interacting pairs. The density of galaxy pairs at these distances 
might by substantially higher than the local galaxy density which 
may produce an excess in the LF at high luminosities. Several au- 
thors find this excess at high luminosities or high masses. For exam- 
ple Lawrence et al. (1999) find a similar excess in their 60 [im LF, 
as do Garcia-Appadoo, Disney & West (in preparation) for their 
HI Mass Function, who find that the higher HI masses are typi- 
cally multiple systems. One can also think of ways our use of a 
global colour-colour relation might have produced a spurious ex- 
cess if, for example, the galaxies at the highest luminosities have 
systematically different colours. This would not, however, explain 
the excess seen in the 60 (xm LF. 

In our earlier work the 850 \im LF derived from the IRS sam- 
ple (DOO) was found to have a slope steeper than —2 at the low 
luminosity end, suggesting that the submillimetre sky should be in- 
finitely bright (a submillimetre 'Olbers' Paradox'). Using the OS 
sample we find the slope of the PSCz-extrapolated 850 (im LF is 
— 1.38, showing that the LF does flatten out at luminosities lower 
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than those probed by the IRS sample, thus solving the submillime- 
tre 'Olbers' Paradox'. 



6 CONCLUSIONS 

Following our previous SCUBA survey of an //MS-selected sample 
of galaxies we have carried out the first systematic survey of the 
local submillimetre Universe free from dust temperature selection 
effects - a submillimetre survey of a sample of 81 galaxies selected 
from the CfA optical redshift survey. We obtained the following 
results: 

(i) We detected 52 out of 81 galaxies at 850 iJ.m and 19 galax- 
ies at 450 ^.m. Many of these galaxies have 850 |j.m emission which 
appears extended with respect to the DSS optical emission, and 
which appears to correspond to very faint optical features. 

(ii) We fitted two-component dust spectral energy distributions 
to the 60, 100, 450 and 850 |j,m flux densities for 18 of the galax- 
ies which were detected at 850 |i^m and at 450 ^.m. We find that the 
IRAS and submillimetre fluxes are well-fitted by a two-component 
dust model with dust emissivity index /3=2. The tight and fairly 
constant ratio of S'45o/S's5o for both the OS galaxies and the IRS 
galaxies is evidence that /3 ~ 2. The temperatures of the warm 
component range from 28 to 59 K; the cold component tempera- 
tures range from 17 to 24 K. 

(iii) We find the ratio of the mass of cold dust to the mass of 
warm dust is much higher for our optically-selected galaxies than 
for our previous work on //?A5-selected galaxies (DEOl), and can 
reach values of ~ 1000. By comparing the results for the IRAS- and 
optically-selected samples we show that there is a population of 
galaxies containing a large proportion of cold dust that is unrepre- 
sented in the IRAS sample. 

(iv) We also fitted single-temperature dust spectral energy dis- 
tributions (to the 60, 100 and 850 |j.m flux densities) for the 41 
galaxies in the OS sample with detections in all 3 wavebands. The 
mean best-fitting temperature for the sample is Td = 31.6 ± 0.6K 
and the mean dust emissivity index is /3 = 1.12 ± 0.05. These val- 
ues are significantly lower than for the IRS sample. The very low 
value of /? is additional evidence that galaxies, across all Hubble 
types, contain a significant amount of cold dust. 

(v) Using our isothermal fits we find a mean dust mass 
Md = (2.34 ±0.36) X 10'' Mq, which is comparable to that 
found for the IRS sample. However, using our two-component fits 
we find a mean dust mass a factor of two higher. 

(vi) We find little change in the properties of dust in galaxies 
along the Hubble sequence, except a marginally significant trend 
for early-type galaxies to be less luminous submillimetre sources 
than late-types. 

(vii) We detect 6 out of 1 1 ellipticals in the sample and find 
them to have dust masses in excess of 10^ Mq. It is possible, how- 
ever, that for some of these galaxies the submillimetre emission 
may be synchrotron emission rather than dust emission. 

(viii) We have derived local submillimetre luminosity and dust 
mass functions, both directly from the optically-selected SLUGS 
sample and by extrapolation from the IRAS PSCz survey, and find 
excellent agreement between the two. By extrapolating the spec- 
tral energy distributions of the IRAS PSCz survey galaxies out to 
850 H-m we have probed a wider range of luminosities than probed 
directly by the SLUGS samples. We find the LPs to be well-fitted 
by Schechter functions except at the highest luminosities. We have 
shown that, whereas the slope of the //fAS-selected LF at low lumi- 
nosities was steeper than —2 (a submillimetre 'Olbers' Paradox'), 



the PSCz-extrapolated LF, as expected, flattens out at the low lumi- 
nosity end and has a slope of —1.38. 

(ix) We find that as a consequence of the omission of a popula- 
tion of 'cold' dusty galaxies from the IRAS sample the LF presented 
in our earlier work (DOO) is too low by a factor of 2, and the DMF 
by a factor of 4. 

In order to further investigate the properties of dust in galax- 
ies follow-up optical imaging (to obtain deeper images than avail- 
able from the DSS) for the whole OS sample detected at 850 H-m is 
needed, in order to make a full comparison of the optical versus 
submillimetre emission. This is important since for many of the 
OS sample galaxies the 850 (im emission appears extended with 
respect to the DSS optical emission. Work on obtaining this data is 
in progress. 
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